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THE GREY TIN WHITE TIN TRANSITION TIN-MERCURY 


SMITH 


ABSTRACT 


The characteristics the transition between white tin (8) and the grey 
modification (@) zone-retined tin containing various amounts mercury have 
been investigated. From dilatometric measurements, shown that the transi- 
tion temperature the reaction: 


alloy) grey tin mercury 


—8.15°+0.1° Above this temperature, mercury will attack grey tin form 
white tin amalgam. The compound previously reported not undergo 
transformation, shown transform below this temperature. Measurements 
are given the rates growth islands grey tin the surfaces number 
tin-mercury alloys —30° These indicate that the rate transforma- 
tion virtually independent the mercury content the alloy. 

The results are discussed and mechanism advanced account for the large 
displacement the transition temperature observed when mercury 
alloyed with tin contrast the small changes brought about other 
additions. 


INTRODUCTION 

transformation tin because mercury will react with white tin for 
considerable interval below the transition temperature. Thus, number 
observations the transformation characteristics this system have been 
made the author during wider studies this transition. 

Earlier, Groen and Burgers (1954) reported the formation compact pieces 
grey tin contrast the fragmentary product usually obtained when 
tetragonal white tin transformed into the diamond cubic grey modification. 
Further investigation revealed that this compact nature had resulted, part, 
from mercurial contamination (Groen 1954; Hall 1955). Groen also reports 
that the hexagonal phase, formed 8-tin when mercury present 
excess 0.2 at.%, unchanged upon transformation the tetragonal 
phase. 

More recently, Ewald and Tufte (1958, 1959) observed that free mercury 
appears upon transformation alloy. This enabled them 
produce single crystals grey tin for semiconductor studies, containing 
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0.001 at.% less mercury deposition from supersaturated solution 
tin mercury temperatures between —20 and —30° They report that 
mercury will suddenly attack these grey tin crystals form white tin solid 
solution the temperature rises above 

Since number the above observations appeared variance with the 
author’s, more thorough study has been made. 


EXPERIMENTAL PROCEDURE AND OBSERVATIONS 


(a) The Transition Temperature 

The transition temperature the system, i.e., that temperature 
which the equilibrium with all its transformation products, was 
investigated follows: 

Alloys were prepared from zone-refined tin* and high-purity 


Alloy 0.1 at.% 
Alloy 6.0 at.% 
Alloy 20.2 at.% 


Alloy was found much harder than whilst was weak 
and brittle, its surface being wet with mercury. 

Three glass dilatometers were used, each containing 5-g specimen with 
xylene the inert fluid (Smith and Raynor 1957). All had 2-mm diameter 
capillary column and, after being sealed, they were calibrated from 
+50° Exposure —27° for few days resulted visually complete 
transformation. Slight crumbling occurred specimen T81 (alloy 1), T82 
showed few large cracks but was little changed shape, and T83 (Hg was 
now composed large sponge-like lumps. Droplets mercury were also seen 
the latter two specimens. 

The volume changes which had accompanied transformation were 
evaluated 10° and are shown Table together with those computed 
for complete transformation pure grey tin and free mercury. Densities for 


TABLE 


Specimen T81 T82 T83 


Mercury content, 6.0 20.2 
density, g/cc 7.31 7.66 8.70 
Specimen weight, 5.065 5.078 5.735 
Measured volume increase, 0.188 0.159 0.147 
Calculated volume increase, 0.188 0.170 0.165 
Measured volume change 24.0 22.2 


original volume 


*Calculated for —10° C. 
+The assigned accuracy of these experiments is 1%. 
$27.1% for pure tin. 


*Vulcan extra pure tin (99.998%) was zone-refined purity 99.9999% (as indi- 
cated transformation rates and analysis). 
Reagent (>99.999% Hg). 
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pure tin 5.75 g/cc for (Thewlis and Davey 1954) and 7.31 g/cc for 
(Lee and Raynor 1954) have been used these calculations. 

Following this, the dilatometers were placed thermostat —4.8° 
the levels the menisci being noted regular intervals. After hours, they 
were returned refrigerator —27° The observed volume contractions 
were interpreted transformation a-tin, T81, 35% T82, and 100% 
T83, that calculated Table This was then repeated other tempera- 
tures, always allowing least hours between successive tests. The results 
are shown Table II. Clearly, mercury begins attack grey tin just below 
—8° and reaction will completion comparatively swiftly tempera- 
tures above —4.8° 


TABLE 


transformed 


Test 
Temp., duration, Order 


Note: All specimens rested at — 27° C for 24 hours between tests. 


establish more closely the temperature which mercury will begin 
react with grey tin, 70-g samples alloys and were carried 
through the above procedure dilatometers with capillary columns, 
but resting period was allowed between tests. Table III contains the results. 


TABLE 


The transition temperature the transformation alloys 


Change meniscus height, 


Test 
Temp., duration, Order T28, T42, 
test 0.1 at.% 6.0 at.% 20.2 at.% 
120 +3.8 +2.1 


Note: Transition temperature is —8.15°+0.1° C. 
*The final positions assumed by the menisci during test 4 were those occupied at this tempera- 
ture at the beginning of these tests. 


Negative signs indicate falling meniscus and transformation, change 
meniscus position being equivalent approximately transforma- 


tion. The transition temperature the reaction: 
alloy a-tin mercury 


—18.5°+0.1°C. 
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The effect mercury content the linear rate phase-boundary displace- 
ment was determined inoculating small sections white tin (approximately 
cm) with grey tin and maintaining them —30° After grey 
tin wart had become established, its boundary positions were noted regular 
intervals enable the rate advance single interface calculated 
for each specimen. These are shown Table IV. 


TABLE 


Rate transformation 


Composition (+10%), mm/hr 
Zone-refined tin (99.9999% Sn) 2.0 
tin 0.1 at.% 0.23 
Zone-refined tin 6.0 0.24 
Zone-refined tin 20.2 0.20 
Vulcan extra pure tin Sn) 1.2 
Vulcan extra pure tin 0.25 at.% 0.24 


seen that the presence mercury zone-refined tin reduces the rate 
transformation considerably, but that this reduction virtually in- 
dependent mercury concentration (for the range studied). Also, the rate 
transformation Vulcan tin, whose principal impurity lead, reduced 
mercury that the other alloys. 


DISCUSSION 


From the tin-mercury phase diagram (Tin Research Institute 1949), alloy 
would appear consist entirely the compound Yet this alloy 
underwent almost complete transformation, thus establishing that 
does transform, with the ejection free mercury. Also, the work 
Section shows conclusively that the temperature which mercury will 
attack grey tin —8.15°+0.1° significantly lower than that reported 
Ewald and Tufte. Absence transformation reported Groen 
(1954) was probably due experiments being carried out above this tempera- 
ture, any mercury rejected during transformation the tetragonal phase 
forming more this compound. (Hexagonal crystals, always form 
from and mercury above (Hansen 1936).) 

interest note that other experiments performed with alloy 
transition temperature 12.9° was observed (Smith (to published) 
value which lies the range expected for dilute tin alloys (Raynor and Smith 
1958). From this and the evidence Tables and III, may concluded 
that any mercury ejected during transformation alloy 
will recombine with some the grey tin form amalgam the tempera- 
ture rises above —8.15° The remainder the grey tin will then begin 
transform the white modification 12.9° The volume changes observed 
tests and Table would further suggest that the average composition 
the new amalgam formed these temperatures corresponds the formula 
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other studies the author (to published), the rate which mercury 
will attack grey tin considerably reduced alloys. droplet 
mercury was placed the surface coherent piece the grey modifica- 
tion wt.% alloy 25° change was observed for 
about hours, after which amalgamation gradually took place. presumed 
that the presence the covalent germanium increased the stability the grey 
tin bonds and reduced the rate attack. Ewald (1954) reported that 0.75% 
germanium will also render grey tin stable +50° 

The results Table are particular interest for, apparently, the 
a-phase advances into alloy, mercury rejected the 
interface. This mercury-enriched boundary layer can only have limited 
width and so, transformation proceeds and further mercury rejected, 
regions free mercury will created behind the interface, probably appearing 
cracks the grey tin matrix. The independence transformation rate upon 
mercury content would suggest that the governing factor the diffusion rate 
tin atoms through the mercury-enriched phase boundary. 

Since the presence lead mercury reduces the rate transforma- 
tion zone-refined tin, the rate the Vulcan extra pure mercury alloy 
would expected lower than the zone-refined alloys. 
The similarity the experimental values would suggest that the inhibiting 
actions lead and mercury are not additive, that the lead being reduced 
the latter. This view supported Groen’s report (1954) that the presence 
mercury 0.03% alloy actually increases the previously very small 
rate transformation. The following mechanism suggested account 
for this. 

Any material present the white tin but insoluble the a-phase will 
rejected the phase boundary and provide barrier its free movement. 
If, however, the inhibiting impurity atom rejected, together with mercury, 
into the mercury-enriched interface, continuing transformation will result 
being transported, diffusion, along the boundary layer until can 
deposited pocket free mercury. Thus, the inhibiting action a-insoluble 
impurities will greatly reduced the presence mercury. 

The expulsion mercury when alloy transforms provides 
considerable clue understanding the transformation tin generally. 
supports the assumption made above that covalent grey tin can dissolve only 
trace amounts solute with metallic bonding. The necessary rejection such 
solute transformation tin alloys would thus tend reduce both 
the thermodynamic instability the white tin when below its transition 
temperature and also the rate conversion upon transformation, both 
which are observed. Also, extremely helpful understanding the 
very small effect that most metallic alloying additions have the transition 
temperature pure tin (Raynor and Smith 1958). Upon initial transfor- 
mation, presumably most the solute rejected and rendered inactive since 
cannot redissolve the white tin phase the temperature raised. 
contrast, rejected mercury, being very mobile even —30° can readily 
recombine the temperature rises above —8.15° form white tin solid 
solution and displace the transition temperature considerably. 
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HEAVY-LIQUID BUBBLE 


ROBERT WILLIAMS? 


ABSTRACT 


Experience with some high-Z liquids fluid mixtures bubble chambers, 
and from other groups, discussed. Measurement procedures are ana- 
lyzed, and optimum methods are presented for angle measurements, deter- 
mination from multiple scattering, and momentum determination from magnetic 
curvature. Magnetic fields will useful even strongly scattering liquids. 
Combined use magnetic and scattering determination not much better than 
one alone. Results are presented detection efficiency for high-energy gamma 
rays, and the severe effects radiation straggling energy measurements 
electron pairs. 


INTRODUCTION 


The bubble chamber was invented only years ago (Glaser 1953) but 
present the principal tool high-energy physics. Propane liquid hydrogen 
has been used most experiments date, but the possibility using liquids 
high atomic number, and therefore good gamma-ray detection properties, 
was early emphasized Glaser, who has constructed xenon chamber 
(Brown, Glaser, and Perl 1956). Other liquids have been investigated the 
Massachusetts Institute Technology (M.I.T.) group (Pless and Williams 
1957; Yamamoto 1957; Jones 1958), the Cal Tech group al. 1957; 
Bugg 1958), and several European groups (Argan 1958; Bertanza al. 
1958; Bullock al. 1958). The primary objective short mean free path 
for gamma-ray materialization. least two experiments the neutral decay 
modes strange particles have already been undertaken with heavy-liquid 
chambers, and more are expected. Figure photograph from run started 
the Cambridge bubble chamber group Harvard, Brandeis, and 
Brown Universities) using propane methyl iodide chamber the Brook- 
haven cosmotron; three pairs and small electron shower 
are 

The present paper examines the properties some practical high-Z liquids, 
and analyzes some the measurements which can made such chamber. 
The high density the liquid greatly reduces the precision momentum 
determination magnetic curvature, without affording multiple-scattering 
measurements comparable with those nuclear emulsion. However, the greater 
stopping power leads more frequent range determinations, and the shortened 
mean free path for nuclear interactions important some experiments. 


received June 1959. 
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experiment will resumed when the cosmotron again operating. xenon exposure 
made Berkeley the Michigan group presently being analyzed. 
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CHARACTERISTICS HEAVY-LIQUID CHAMBERS 


Nearly any liquid which stable temperatures approaching its critical 
temperature would function bubble-chamber liquid; problems tem- 
perature and pressure, corrosiveness, expense, and toxicity limit the choice. 
The theory bubble-chamber operation (Pless and Plano 1956; Seitz 1958) 
indicates that surface tension the controlling property: vapor pressure 
operating conditions not too high, the surface tension should low. 
The surface tension high-Z liquid can lowered admixture 
lighter liquid, even dissolved gas. The properties mixed-fluid chambers 
are not strikingly different from those pure-liquid chambers; they were 
first investigated Argan and Gigli (1956); the application high-Z 
chambers was made Pless and Yamamoto (Yamamoto 1957; Pless and 
Williams 1957), who showed that propane and iodide yield con- 
venient mixture. 

Table gives some characteristics high-Z bubble-chamber liquids which 
have been tested various groups. Xenon and present special difficulties 
procurement handling; the other compounds are articles commerce. 


TABLE 
Properties some bubble-chamber liquids 

Radiation length, Xo, 8.6 13.7 46.5 
Radiation length, 3.7 3.8 10.5 9.5 17.6 105 
100-Mev gamma mean free 


The characteristic length for radiative energy loss extreme relativistic 
particle the radiation length (defined, for example, Rossi 1952); 
also the natural unit for related processes such pair production and Coulomb 
scattering. list the radiation length (including Born approximation and 
atomic-electron corrections) for the various liquids, and also the effective 
mean free path for pair production 100-Mev gamma rays. The latter 
includes the pairs created the field atomic electrons (one the three 
electrons usually very low energy compared with the other two), but 
excludes any pairs for which one electron has too low energy (<4 Mev) 
make recognizable pair. found the discussion Ashkin and Bethe 
(1953) convenient for pair-production cross sections. 

The presence hydrogen the mixture proves very desirable, 
despite the fact that one cannot make good production cross-section measure- 
ments. One finds that some the events (for example, associated production 
and particle certain angular ranges) have overdetermined kine- 
matics that they can definitely identified events, and can therefore 


PLATE 


3 


Fic. Negative pions incident from the left methyl iodide propane 
mixture. Two gamma rays from the end track and one from track create pairs the 
chamber. small electron shower also visible. 


j 
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used determine beam energy and measurement accuracy. thorough 
discussion the separation hydrogen events from carbon events propane 
will found the work the Columbia group (Eisler 1958). have 
used their results the reaction estimate the number 
such associated-production events heavy liquids, and the fraction 
occurring free protons, follows: find the cross section for 
associated production carbon four times that free proton: 
ata energy about Bev. For heavier nuclei assume 
that associated production occurs the same ratio carbon the measured 
inelastic cross sections for high-energy pions (Cronin 1957), that is, that 
the fraction inelastic events yielding strange particles when pions strike 
complex nuclei independent This would not true proton beam, 
but probably fair approximation for pions. Table gives the fraction 
associated-production events occurring free for the various liquids. 


MEASUREMENT ACCURACY 


The mechanism bubble formation bubble chamber very probably 
associated with delta rays low energy (order kev), most which have 
range considerably less than micron. The bubbles grow much too rapidly 
undergo noticeable diffusion, that, the absence turbulence 
statistical asymmetry bubble growth, the center bubble should corre- 
spond the path the particle which seeded better than micron. 
The accuracy track location, Glaser has pointed out (Brown 1956), 
may principle comparable with that photographic plate rather than 
that cloud chamber, although the number locations per centimeter 
only 50, cloud chamber. practice one forced drastic com- 
promises the optics bubble-track recording: the bubbles must grow 
perhaps 100 microns order visible, and small numerical aperture 
must used have manageable photography and reasonable depth 
field. The result that the standard error the location one point 
track may the order 0.005 0.01 cm. The accuracy most critical 
the application multiple scattering momentum determinations; 
multiple-scattering experiment with small chamber and good optics Kirillov 
al. (1958) seem have had standard error (estimated from equation (1) 
below) only 0.002 0.003 cm. 


Angle Determinations 

Most the relevant measurements (kinematic reconstruction, momentum 
determination from magnetic curvature, from multiple scattering) depend 
determining the angle portion the track. The problem, one 
plane, the least-squares fit straight line length number 
measurements, per centimeter, each which has r.m.s. error (We 
use throughout for absolute errors, for relative errors.) This gives angle 
error 


can greater than the bubble density, course, and practice will 
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often less. the error the direction perpendicular the photographic 
plane, and its corresponding dip-angle error, will greater geometrical 
factor depending the stereo system. For two cameras separated distance 
objective distance the factor, assuming independence the errors 
bubble location the two images, the largest usable 
angle the lens one often has D/u, that the factor 
typically about 15°). 

The optimum measurement length do, the absence distortions, will 
determined minimizing the measurement error, above, combined with 
the error due multiple scattering. convenient express multiple- 
scattering angles terms the formula 


where the track length the radiation length, the momentum- 
velocity product Mev, and which varies slowly with 
and with the the material, and course depends numerically what 
angle one measuring. (We assume singly charged particles.) 

Table gives values for the mean projected angle between true 
tangents the ends the track. want instead the projected r.m.s. angle 


TABLE 


Projected mean Coulomb scattering angle radians between tan- 
gents track length material radiation length and 
atomic number for particle velocity Table gives 
the formula 


Mev 


corresponding the least-squares fit straight line the track. Scott 
(1952) gives, for Gaussian statistics, the factors .765; from 
Table find typical value for short fast tracks heavy 
liquid; for our analysis can assume and the 
multiple-scattering error Minimizing the sum the 


squares the errors yields 


the corresponding minimum angle error 


where the units are centimeters, million electron volts, and radians. 
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noteworthy that the dependence the error not strong, 
that the heavy-liquid chambers will only factor two three worse 
than propane chambers this respect. The Columbia group (Eisler al. 
1958) reports average angle errors about .004 radian propane (which 
implies, for 400 Mev, location error .01 cm). When magnetic 
field present there will correction the apparent direction; this 
discussed Section 


Multiple Scattering 

The possibility using multiple-scattering measurements determine 
the bubble-chamber tracks, following the well-known techniques 
nuclear emulsions, has been explored several laboratories; Teem and co- 
workers (Alyea 1957) have discussed results from very small chamber, 
while Kirillov and co-workers (1958) report measurements propane. The 
crucial point here the measurement accuracy, since for arbitrarily good 
accuracy one can obtain estimate (p8c) limited only the total number 
bubbles. convenient consider first the theoretical expression for 
mean-scattering angle, then the actual measurement, and defer the 
next section discussion the problem when magnetic field present. 

Multiple scattering discussed detail Scott (1952); the most con- 
venient statistic the mean absolute projected angle (the actual track 
projected onto plane containing the initial direction the particle, and 
the angle between the initial direction and the projection the final 
direction considered). Within the restrictions small-angle approximation 
the full distribution (not Gaussian approximation) derived, and the mean 
angle calculated under the assumption that any single scattering more 
than four times the average scattering for the whole track will rejected— 
Scott calls this the cutoff case. For very high momenta and long tracks the 
finite size the nucleus causes decrease scattering not calculated 
Scott. Cooper and Rainwater (1955) have given numerical results for the 
scattering distribution including nuclear-size effects, but bubble-chamber 
scattering measurements will not ordinarily extend such high-momentum, 
long-track cases. Rejection large-angle scatterings (the cutoff procedure 
described below) greatly reduces the error caused inaccurate treatment 
the tail the distribution. 

Scott’s numerical results (or, equivalently, those Moliere which agree 
very closely) for the projected mean angular deflection, particle 
traversing distance are presented table values (equation (2)) 
Table which has been taken from the larger table Barkas and Rosen- 
feld (1958). These values are for the non-cutoff case. 

estimate (p8c) from multiple-scattering measurements one would like 
the true direction every point along the track. One has available the 
positions certain points, within certain error (the assumed 
Gaussian with standard deviation and one must integrate out the noise 
(and subtract out what not integrated out) obtain the average changes 
direction, from second differences some comparable operation. The 
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following method, which must fairly close optimum, frequently used 
nuclear emulsion work‘ and will serve for analysis. 

One divides the track length into equal segments, and makes 
least-squares fit straight line each track segment. The angles 
between adjacent lines then give the observed mean scattering 1/n 
Any angle which more than four times the mean angle 
then excluded from the data and the mean, computed. The 
remaining contribution scattering will include the error given 
equation (1); may also include appreciable contribution from the measuring 
engine and from turbulence. The mean noise scattering, for given cell 
length can determined from measurements high-momenta 
The best estimate for the true mean angle scattering 
will 


Table gives the mean change direction path length how- 
ever, the (cutoff) mean change direction between two best-fit lines, each 
length Scott presents results for related case, the cutoff mean angle 
between two chords (Scott 1952, equation (8f); note that 2d) and states 
that the effect least-squares fit can included multiplying his results 
for 1.05. terms Scott’s function the scattering 


and the multiple-scattering estimate 


compare with Table compute representative for bubble-chamber 
from Table II, 10.4 Mev. 

find optimum value is, appropriate noise-integrating 
length—we use equation (1) for single segment r.m.s. angle get mean 
difference noise angle given 


say mo, exist which minimizes the relative error and therefore 
The nature the noise distribution and the accuracy with which 
known will influence one’s choice the best ratio for photoplate 


‘In nuclear emulsion work, however, chord-fitting may done eye. Bubble-chamber 
tracks are usually measured semiautomatic co-ordinate recorder and the line segments 
would fitted computer. 

(1950) gives method for subtracting the noise using the internal statistics 
the track; essentially, differs from zero only because correlations introduced 
noise. The method the text should superior good estimate available. 


4 a 
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work the problem discussed Dilworth al. (1950) and Voyvodic 
(1954). assume that accurately known and its fluctuations are 
Gaussian, then from (5), (7), and (8) can find the optimum ratio 
tice would suggest greater cell length, lessen the contribution but 
for bubble-chamber accuracy this should serve well enough. From (7) and 
(8), taking 8.5 Mev have 


(0.8 

Loe, pBe 
which essentially the same (3). Even crude guess for will serve 
estimate 

The relative standard deviation the mean-scattering angle can esti- 
mated from the statistical discussion Scott (1952), who finds about 
there noise; with our noise ratio this increased (1+1/7) give 

0.97 
(10) 


ever, the distribution function not Gaussian, and any case its 
inverse, for typical values will decidedly asymmetric. fairly 
optimistic example consider track 20cm long with 300 Mev, 
increase track length appears difficult visualize dramatic improve- 

practice corrections must made for the dip the track (since 
usually more accurate and convenient project onto plane parallel the 
film plane rather than plane containing the beginning the track), and 
for energy loss along the track. Although scattering increases the particle 
loses energy, very little improvement statistical accuracy obtained 
progressive shortening the cell length. make energy-loss correction 
one must either estimate assign mass the particle. Using the 
mation (Annis a/. 1953) const. one finds first approximation 
that the estimate equation (7) should raised give 
the the beginning the track. (Ro the range corresponding the 
initial energy.) 

Subtraction magnetic curvature considered Section Experimental 
tests the multiple-scattering theory nuclear emulsion, 
references, can found the paper and Eisenberg (1958). 

case particular interest the identification particle which stops 
the chamber, estimating its mass from multiple scattering. The problem 
has been solved Olbert and collaborators (Annis 1953), who give the 
distribution function for variable related the scattering angles and 
residual range, measured points along the track. For single track 
reasonable (20-30 cm) length the separation between pions and protons, for 
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example, not very good: the small-scattering part the pion distribution 
(it chi-squared curve with degrees freedom) overlaps the proton 
distribution considerably, that tracks with little scattering are ambiguous. 
Widely scattered tracks, however, can identified pions with good confi- 
dence. The appropriate cell length can estimated applying 
equation (3) the beginning the track; little gained decreasing the 
cell length toward the end the track. The mass (half-width 
point) about Am/m 1.6/+/n. One sees from equation (3) that the 
high the liquid not great advantage multiple-scattering measure- 
ments, since the standard deviation varies only 

homogeneous group stopping particles can lead fairly good estimate 
their mass the Olbert method; results for propane chamber are given 
Kirillov al. (1958). 


Magnetic Curvature 

magnetic field modest magnitude (15-20 kilogauss, say) can lead 
useful momentum measurements chamber even quite short radiation 
length and course one can visualize much larger fields with pulsed 
magnets. The general problem extracting the maximum information from 
combined magnetic and multiple-scattering measurements has been con- 
sidered number authors: the most explicit and useful result the 
work Kim (1958). However, the results are disappointing the sense that 
one finds little improvement over the usual practice choosing only the better 
the two estimates. 

Kim proves the validity making the individual magnetic and multiple- 
scattering estimates and taking the weighted average. This leads 
improvement the error, factor which always less than 
V2. numerical example which each estimate separately has minimum 
error 10%, Kim finds that the combined error 8%. 

shall therefore confine this discussion the individual estimates. 

Correction Multiple Scattering 

The magnetic field causes non-zero mean angular deflection which must 
estimated finding the net angular deflection 


n 
= 


the scattering angle then 


and all goes through before. 


Magnetic Estimate Momentum 
The momentum component particle, transverse uniform magnetic 
field related the angle through which its projected path would 


turn the absence scattering 
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where Mev, gauss, and the length projected path centi- 
meters. well known that the best estimate for track curvature, the 
presence fixed amount multiple scattering and noise, obtained from 
the observed net deflection the track (see, for example, Scott 1949; 
Moyal 1950). Random scattering, other words, automatically averages 
itself out the most efficient way and amount inspection the shape 
the curve can improve this. One normally estimates sub- 
stituting for (11), with standard deviation given the root- 
mean-square multiple-scattering angle, estimated (using Gaussian 
statistics) from (2) from However, Kim (1958) shows that 
the estimate smaller than when scattering 
prominent; gives graphical solution but for practical cases his approxi- 
mate relation adequate: 


Arms 2) 
(12) (Omag est 1— | 
mag 


where the relative standard error the magnetic estimate, given from 
(2) with Mev, and (11), 


(13) 


Here gauss, the true track length centimeters, and the dip 
angle. Exclusion large single scatterings will reduce 15%. 

The correct estimate for the imitial direction track must include the 
effect the mean curvature the track, this applies the 
direction measurements Section well the estimation the 
total angle least-squares fit parabola known curvature toa 
segment track length shows that the slope given the straight-line 
least-squares formula must corrected The corrected estimate 


for the total angle therefore 


(14) Qeorr 


and this should used (12), iterated necessary. 
The ratio magnetic multiple-scattering error can used indicate 
which method preferable; for non-dipping tracks 


the usual units, and independent track length. With 0.005 cm, 
and our previously cited example 300 Mev, 

Under our assumptions turbulence and excellent optics the uncertainty 
magnetic determination will usually have contribution from measure- 


ment errors which less than half the error due scattering, rm. 
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From (4) and (13) the ratio may written 


Energy loss along the track will often important correction. one 
assumes power law relation between momentum and range, and assigns 
mass, one can find the initial momentum from the observed total angle; 
first approximation (11) gives the momentum the center the track, and 
(13) correct refers the center the track. 

The latter part track which near the end its range does not add 
much information momentum determination. practical cutoff point 
can obtained from equation (3) requiring that the fractional energy 
change distance not exceed some reasonable value. 


Sign the Charge 

When large one can still ask the sense the net deflection gives 
reliable indication the particle’s charge; this problem was solved Moyal 
(1950). Both the total deflection and mean-scattering angle must 
measured, and estimate must somehow made, independently 
(one might estimate from and guess from bubble density, 
mass assignment, one might have limits from external con- 
siderations). Then the confidence level the sign indicated (that is, 
the fraction such observations which would give the right answer there 
were equal priori probability for either sign) 


Measurements Electrons 

obtain information high-energy gamma rays one must look 
electron pairs; the opening angle too loosely correlated with gamma-ray 
energy, and usually too small, much use, one has estimate the 
electron energies. The accuracy this measurement, high-Z liquid, 
severely limited the uncontrollable effects radiation straggling. The 
large fluctuations energy loss radiation are familiar phenomenon 
electron measurements; illustrate point out that high-energy electron 
traversing one-half radiation length will suffer average energy loss 
one-half, but most probable loss zero; the straggling function rises mono- 
tonically infinity zero energy loss. 

There are two types measurements might consider. The first the 
determination parameter known distribution from observations 
the energies sample electrons. The procedure straightforward: one 
“folds” the distribution with the appropriate radiation straggling function, 
and determines the parameter from the folded distribution maximum- 
likelihood equivalent methods. This method has had number appli- 
cations; described, for example, Vilain and Williams (1954) with 
reference the electron spectrum from muon decay, and was applied 
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bubble-chamber measurement Pless al. (1957); shall not discuss 
further. The second, and usually more pertinent, problem the estimation 
the energy individual electron, which discuss this section. 

The probability, that electron energy have 
energy after radiation lengths has been given Eyges 
(1949) various approximations, which the one appropriate energy 


are neglecting collision loss. What want inverse problem: given 
observed energy (from magnetic curvature, say), with what confidence 
can say that this electron’s original energy lay within some specified 
range; that is, what fraction out all electrons that give would 
come from that range the experiment were repeated indefinitely? Evidently 
this may depend the say all possible 
electrons which might observed the experiment. The usual situation 
that distribution functions for errors, collision loss, etc., are well behaved 
that the answer our question essentially independent the form the 
priori distribution shall show that this mot the case for radiation 
straggling, taking the physically form for the 
uniform distribution per logarithmic energy interval, 


could measure the energy the electron after path length 
could use, our inverse problem, the fraction all electrons energy 
which come from the range E+AE, which clearly 


E+AE 
E E 


However, must use the entire path for our measurement, and from the 
total track deflection, which the integral the reciprocal the energy, 
obtain the observed energy make progress assume that the 
entire energy loss occurs the midpoint the whence 


(19) 2/Eops = 1/Eot+ W/E. 
Now can calculate the fraction all electrons energy which 


Eopst 


Eovs 


obs 


where must eliminated (19). Using (17), (18), and (19), can 
express terms the incomplete gamma function, tabulated Pearson 


can see that more accurate treatment would give the same results the lowest 
order 


— 
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where have put This not even approximately Gaussian; 
makes 0.68, say That is, there 68% chance that lay within 
Which the confidence level one standard deviation the 
Gaussian case. 

gives the confidence level 95%, which would two standard deviations 
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FOR 


CONFIDENCE LIMIT, 


PATH LENGTH ,t 


Fic. Radiation straggling effect. The confidence limit plotted against 
path length radiation lengths, for 68% and 95% confidence level. Dashed lines show 
the standard deviation magnetic curvature measurement for two values Blo. 


the Gaussian case; course this case much greater than pick 
the path length which would give the optimum estimate must 

compare with the error magnetic curvature from (13). Figure 


Large they are, these errors are deceptive because they contain contri- 
bution (though less than half) from the radiation fluctuation. ask for 
the 95% confidence limit, with the better case cos find 


05+ SS + + 
0.4} + + + 
| 
! Trad , 68% | 
| | 
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different and much shorter optimum, 0.14, 1.5. 
measurement multiple-scattering methods presumably would give the 
same result. 

show the degree which the entire calculation dependent the form 
the priori distribution have calculated for cutoff distribu- 
value given Fig. i.e., 0.23 0.20. 

The above errors will course multiplied factor for gamma- 
ray energy and factor which can large for energy 
(four electrons). clear that compromise would necessary 
one wanted measure gamma-ray energies better than 50%. the 
next section the effect gamma detection efficiency considered. 


Gamma-Ray Detection Efficiency 

The strongest demand for short radiation length (high average set 
gamma-ray materialization probabilities. Table gives the effective pair- 
production mean free path for 100-Mev gamma rays, about 1.8/) (for 
400 Mev would about here present rough calculation the 
efficiency some typical observations, assuming constant value for 

assume spherical observation volume, radius with the gamma-ray 
source distributed uniformly throughout it, and isotropic, non-correlated 
space distribution for the gamma rays. (The latter partly realized 
mesons and hyperons reasonable energy which decay via neutral pions, 
for example.) Let the average probability detecting one gamma 
ray which originates want the average over the sphere. 
avoid the unpleasant integral will approximate algebraic 
expression which has the correct value and slope and and 
which approaches the correct limit (unity) plausible way approaches 
infinity. turns out that the edge the sphere can easily calculated 
and 

center course 
with zero slope that point. The basic approximation consists representing 
P(r) for small quadrant ellipse, then replacing the portion from 


the uniform value for arbitrary The result 
integrating this over the sphere gives the average probability see 


single gamma ray: 
(24) 


two gamma rays are emitted (as, for example, the decay neutral 
pion) the probability detect them both (by our independence assumption) 
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and this can averaged over the sphere give similar (24) 
but more complicated. felt that the approximations are good about 5%. 

quantity, also plotted; this the probability seeing 
least one gamma from two-gamma event. For comparison the simple 
central probability also shown. 


OF MATERIALIZATION 


PROBABILITY 


= 


0.24 6 8 | 121416 18 2 222426 28 3 32 34 36 
RADIUS IN MEAN FREE PATHS, R/A 


Detection probability for gamma rays isotropically released throughout spherical 
volume. refers one gamma ray, both gamma rays pair, and least 
one pair. For comparison, the dashed line gives the probability the center the 
sphere. 


good detection probability—requisite for many problems hyperon and 
K-meson decay—evidently demands that the radiation length should 
considerably smaller than the the chamber. Section saw 
that this precludes accurate measurements gamma-ray energies, unless 
one goes very large chamber very high magnetic fields. 


SUMMARY 


The heavy-liquid bubble chamber can provide good (~50%) mean detection 
efficiencies for high-energy gamma rays and neutral pions over large volume. 
The large mass material has other advantages which depend the experi- 
ment question: charged mesons and baryons considerable energy (typically 
~100 Mev) will stop the liquid, with consequent accurate energy 
determination, and observation decay products; nuclear interaction proba- 
bilities are also relatively high. The enhanced multiple scattering should 
permit estimates (p8c) the tedious method familiar from nuclear emul- 
sions, but with lower accuracy than emulsion work. magnetic field will 
give momentum estimates that, within typical range parameters, can 
either better worse than multiple scattering; any case they will 


| 
= 
og | . + + + > + 
H | 
| ic P, | 
| | 
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easier, more reliable, and endowed with the sign the charge. Multiple- 
scattering effects not seriously interfere with geometrical (angle and 
position) measurements. 

Measurements electron energy present separate problem, due 
fluctuations energy loss radiation. detailed consideration this 
problem shows that any reasonably heavy liquid effectively prevents energy 
measurements better than ~50% with present techniques. 
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ULTRASONIC ATTENUATION SUPERCONDUCTING AND 
NORMAL 


ABSTRACT 


Ultrasonic attenuation superconducting and normal mercury has been 
studied with 20, 34, and Mc/s longitudinal waves. The difference attenua- 
tion between the superconducting and normal states below the transition tem- 
perature greatly increased the higher frequencies. The attenuation the 
normal state the same for both longitudinal and transverse critical magnetic 
fields and independent field strength between the critical value and 880 
oersteds. The results disagree some respects with those MacKinnon and 
Myers. The experimental limitations the method critical test the 
theory superconductivity are discussed. 


INTRODUCTION 


Ultrasonic attenuation superconducting lead, tin, indium, and aluminum 
1954; MacKinnon 1955, 1957; Mason and 1956; Morse, 
Tamarkin, and Bohm 1956; Chopra and Hutchison 1957) has been found 
diminish rapidly when the temperature lowered below the transition tem- 
perature. Recently, MacKinnon and Myers (1958) have 
observation the case superconducting mercury. mercury readily 
available very pure state, one expects find marked difference elec- 
tronic absorption ultrasound the normal and superconducting states. 
These authors have observed quite different changes attenuation with 
temperature two specimens for one which the purity specified 
99.999%. The observations were carried out Mc/s. this frequency 
significant diffraction loss occurs addition the two unknown seal losses 
arising from the mercury copper and copper quartz end junctions 
employed. These losses are probably significantly greater than the attenuation 
the superconducting mercury. have eliminated the copper disks and one 
seal and have been able work the higher frequencies 20, 34, and Mc/s 
thus reducing the diffraction losses. MacKinnon and Myers have also observed 
unusual increase attenuation the lower temperature with transverse 
magnetic field. The attenuation the normal state, obtained applying 
the critical transverse magnetic field, has been found much higher than 
that obtained applying the critical longitudinal field. These results are 
unusual and need verification. 


EXPERIMENTAL ARRANGEMENT 
Measurements ultrasonic attenuation were made the conventional 
pulse technique. Pulses (about 1.5-microsecond duration) the required 
frequency oscillation were fed rate about 1000 per second 1/2-in. 
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X-cut quartz crystal which converts them into pulses longitudinal waves 
the same frequency. The quartz crystal the end the cylindrical specimen 
acted transmitter the initial pulse well the receiver the successive 
echoes. The pulse and its echoes were amplified preamplifier and wide 
band amplifier Mc/s passband) and exhibited oscilloscope 
(Tektronix 514D). The pulses were photographed their amplitudes measured 
means scale and cathetometer. The ultrasonic attenuation was 
calculated from the rate decay the successive echoes. 

The cylindrical specimens were prepared freezing mercury cylindrical 
bakelite container inner diameter mm. Four specimens length varying 
from were used different runs. The mercury used for these 
specimens was supplied the Fisher Scientific Company and contained less 
than 0.0005% gold and silver and 0.0001% base metals. After the quartz 
crystal was wetted with distilled water, was placed the liquid mercury 
specimen. The specimen container floated small pool mercury. was 
fixed position first freezing the pool. The mercury specimen was then 
slowly cooled with its top and bottom surfaces parallel and horizontal. The 
thin layer ice sealed the quartz crystal the mercury specimen. The ice 
seal has been found quite satisfactory low temperatures with suf- 
ficiently low seal loss enable make measurements several echoes 
Mc/s. 

The cryostat used was conventional design. Temperatures down 
1.2° were obtained pumping and maintained Cartesian manostat. 
The temperature was determined measuring the vapor pressure over the 
liquid helium bath. Temperatures 4.6° were obtained pressurizing 
the helium Dewar. 

Magnetic fields were obtained two liquid air cooled Helmholtz coils; one 
was used for longitudinal fields and the other produced field perpendicular 
the axis the Dewar and the specimen. The highest field available 
during this work was 880 oersteds. 


RESULTS 


Figure shows plot the attenuation Mc/s longitudinal waves 
pure mercury function temperature. The ‘normal state’ curves were 
obtained applying the critical magnetic field along the direction the 
specimen (longitudinal) and perpendicular the direction the specimen 
(transverse). The critical fields used were taken from the known threshold 
curve for mercury. remove the uncertainties introduced the frozen-in 
flux, the specimen was warmed above the transition temperature after each 
pair superconducting and normal measurements. 

Figure shows plot 20, 34, and Mc/s longitudinal 
waves pure mercury function temperature. The dotted curves are 
those the normal state obtained applying the critical longitudinal 
magnetic field. 

study the dependence ultrasonic attenuation the transverse 
magnetic field, the attenuation various temperatures plotted against 
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Fic. Temperature dependence ultrasonic attenuation Mc/s longitudinal waves 
superconducting and normal states pure mercury. 

Fic. Temperature dependence ultrasonic attenuation 20, 34, and Mc/s longi- 
tudinal waves superconducting and normal states pure mercury. 


transverse magnetic field shown Fig. These plots are quite similar 
those obtained applying longitudinal magnetic field. The latter are not 
shown here. The increase attenuation with transverse field reported 
MacKinnon and Myers not observed. 
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DISCUSSION 


The temperature dependence ultrasonic attenuation the superconduct- 
ing state similar that observed other superconductors. The dependence 
frequency the difference attenuation the normal and superconduct- 
ing states quite marked. Our results show that for frequencies below Mc/s, 
very small change attenuation will observed. 

the normal state, the attenuation increases rapidly the temperature 
lowered. The increase the same for both longitudinal and transverse fields. 
Furthermore, the attenuation independent the applied magnetic field 
(both longitudinal and transverse) long greater than the critical field 
value 880 oersteds. Measurements higher magnetic fields are 
desirable. 

Recently, attempts have been made (Morse and Bohm 1957) use ultra- 
sonic attenuation measurements verify the 
(BCS) theory superconductivity (1957). According this theory, the ratio 
attenuation the superconducting state that the normal state 
given by: 


where the energy gap the Fermi surface and has been calculated 
Bardeen. Apart from the question the validity the above-mentioned 
equation, there little doubt that the ultrasonic attenuation measurements 
date are too crude able distinguish between this relation and 
power law. The superconducting attenuation (a,) measured includes losses due 
the seal (between crystal and specimen), diffraction, non-parallelism the 
ends the specimen, etc. Over the temperature ranges studied, these losses are 
assumed constant. Another uncertain factor the normal attenuation 
a,. two observers have obtained the same dependence with tempera- 
ture. fact, whereas the measurements single crystals lead 
(1954) and tin Mason and Bémmel (1956) show marked dependence 
temperature, which somewhat similar that observed here Fig. 
for pure mercury, those Morse and Bohm (1958) single crystal tin and 
pure polycrystalline indium show little dependence temperature. This 
probably because the measurements are quite sensitive the purity 
the specimen. The values will also affected the presence frozen-in 
magnetic flux determination follows immediately measurement 
With the normal state produced magnetic field. The presence this 
error difficult assess from the literature since few authors have reported 
the details their experimental procedures. 

The results Fig. indicate the marked rise with decreasing tem- 
perature high degree purity the specimen. The values also seem 
approach ‘background’ value expected the seal losses, etc. are constant. 
interesting note, perhaps not too significant for the reasons given 
above, that the theoretical curve calculated the basis the BCS theory 
predicts less rapid fall attenuation, with temperature below than 


— 
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here experimentally found. The comparison made Fig. Morse and 
Bohm (1957) have observed similar qualitative disagreement with BCS 
theory the case superconducting tin and indium. 


BCS THEORY 


EXPERIMENTAL DATA 
FROM FIGURE | 


Fic. Theoretical and experimental variation against 
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FIRST ADSORBED LAYER NITROGEN PYREX 77.4° 


ABSTRACT 


measure the nitrogen adsorption isotherm pyrex (Corning 7740) 
77.4° has been made for pressures above the adsorbed layer ranging from 
Hg. Ultrahigh-vacuum techniques were used. Geometric 
adsorbing areas were and Coverages ranged approximately from 
0.001 0.3 monolayer. Full equilibrium was probably not achieved but the 
data could represented the Dubinin—Radushkevich equation, 


for from where the amount adsorbed. 


INTRODUCTION 


Experimental data the physical adsorption nitrogen below pressures 
and temperatures near the normal boiling point are scarce, 
even though ionization gauges readily operate below this range (Dushman 
1949; Alpert 1953). Further isotherm data for pressures below 
appeared desirable because: 

(a) Surface coverages are less than one monolayer this range and the 
experiments should provide decisive tests for basic theories physical adsorp- 
tion (de Boer 1956). 

(b) Low coverage data has been related the heterogeneity solid sur- 
faces (Halsey 1952). 

(c) Recent theories chemisorption include essential part the van 
der Waal’s interaction between gas and surface (Kisliuk 1957). 

(d) Calculations the thermodynamic functions adsorption (Drain 1954) 
may readily made isotherms are taken different temperatures. 

(e) detailed knowledge the physical adsorption nitrogen desirable 
because its wide use the measurement surface areas (Joy 1953). 

(f) The design vacuum systems containing refrigerated traps may 
facilitated such low pressure data. 

was observed Hobson and Redhead (1958) that physical adsorption 
surface areas the order cm? was readily detectable systems 
volume about liter very low pressures (approximately Hg). This 
observation was extended give measurement the helium adsorption 
isotherm very low pressures 4.2° pyrex (Hobson 1959). The present 
work applies the same exploratory method the physical adsorption 
nitrogen. 

With nitrogen gauge (Bayard and Alpert 1950) was used 
the pressure-sensitive element because its low pumping speed was useful. 
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This gauge pumps measures (Young 1956; Bills and Carleton 1958) and 
may re-emit some the gas which has pumped (Varnerin and Carmichael 
1957). However, these gauge processes had not been examined the above 
authors the pressure range used the present measurements, and for 
unambiguous interpretation some the adsorption data, least semi- 
quantitative knowledge these processes was required. The gauge processes 
have been examined concurrently with the adsorption measurements and 
paper these processes will submitted for publication shortly. 


II. APPARATUS 


The apparatus shown Fig. almost identical with that reported 
previously for helium adsorption (Hobson 1959). The envelope was pyrex 


MAIN BAKING REGION 


BAG BAG 
VALVE 
MMERSIBLE B 2) C AZ) 
FINGER 
MOVABLE 
DIFFUSION 
PUMP 


Fic. Apparatus. 


glass (Corning 7740). The valves were ultrahigh-vacuum valves (Alpert 1953; 
Bills and Allen 1955). The gauges (B.A.G. Nos. and were 
standard design and were operated with collector ground potential (i.e., 
zero); filament, +45 volts; grid, +300 volts; and grid current variable be- 
tween milliamperes and microamperes. The ion current was 
continuously recorded currents amperes. The main adsorption 
volume, (938 was the left valve and contained magnetron 
pressure gauge not used the measurements, gauge (B.A.G.- 
which was the major measuring element the experiment, and finger 
which could immersed liquid nitrogen 77.4° Use double Dewar 
reduced the loss rate the inner Dewar low value. The level 
liquid the inner Dewar was visible all times and could readily adjusted 
fixed mark the finger. Valve controlled nitrogen flow into Valve 
controlled flow from the gas source spectroscopically pure nitrogen 
the central volume, V2, where pressure, P2, was measured B.A.G.-2. 
could pumped through valve with mercury diffusion pump suitably 
trapped with liquid nitrogen. 

The system was prepared for measurement baking the region shown 
Fig. for hours 450 400° with all valves open and with the mercury 
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diffusion pump running. Before measurements were begun pressures 
were obtained with valve closed and B.A.G.-1 operating 
milliamperes. This background pressure mainly the result the gauge 
pumping against the helium diffusing through the pyrex envelope (Rogers, 
was obtained with valve open prior breaking the seal the gas source. 


EXPERIMENTAL PROCEDURE 


typical adsorption isotherm run was performed follows. After B.A.G.-1 
was outgassed with valves and open, was closed and equilibrium 
obtained about minutes. The finger was then immersed fixed mark 
with the grid current milliamperes keep the helium 
background pressure low. This immersion caused pressure drop about 
was next reduced microamperes (or microamperes) reduce 
the pumping speed the gauge during adsorption measurements, and was 
left this value. Upon reduction the pressure began rise slowly 
because the pumping speed the gauge helium was also reduced. The rise 
was recorded until its characteristics were established (after minutes) 
and the pressure reading any future time due helium diffusion could 
predicted. After valve was closed steady high nitrogen pressure was 
Valve was then cracked open and nitrogen allowed bleed into until 
pressure about Hg, about five times the estimated back- 
ground pressure, was established. Valve was then closed and allowed 
stabilize for minutes. The pressure excess estimated background 
was taken the nitrogen pressure above the layer adsorbed the inside 
surface the finger. The Dewar was then lowered, permitting the adsorbed 
nitrogen desorb. The pressure reached peak, from whose magnitude the 
quantity adsorbed was calculated, and then the gas was pumped 
down until the characteristic pumping time order minutes 
could established. There followed set operations shown Fig. This 
set was repeated increasing pressures pressures Hg. With 
reference Fig. may seen that the pressure fall when the finger was 
immersed (characteristic time was much faster than when the gauge was 
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pumping alone; that is, This relation held throughout the pressure 
range investigated and indicates that the slowly pumping gauge represents 
minor fraction the total pumping speed the system the presence 
liquid nitrogen. The variations observed will discussed the paper 
gauge processes. 

Figure was analyzed with the following relations: 


where the gauge ion current, 
the gauge electron current, 
the gauge constant (mm for the gauges used. 
Over most the pressure range covered the pressure above the adsorbed 
molecular formula, 


For from Hg, eq. (2) was corrected with the formula 
Liang (1951). 
Once valve was closed was assumed that 


where the number molecules physically adsorbed, 
the number molecules gaseous phase, 
the number molecules ionically pumped the gauge since 
closure 
Now 
(4) 


sary minor correction was applied eq. (4) because small portion 
was not room 


G 


Combining eqs. (1) and (5), 


where the ionic charge coulombs obtained integrating the 
observed current following closure 

When liquid nitrogen removed from the finger reduced zero. Since 
and are given eqs. (4) and (6), respectively, eq. (3) may then 
found. Once known, can found any time combining (3), 
(4), and (6). The adsorption isotherm concerned with the simultaneous 
value and The values given Section were obtained just before 
liquid nitrogen was removed from the Dewar. 
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The adsorbing area was taken the simple geometric area the inside 
surface the finger (inside diameter, 2.2 cm; outside diameter, 2.5 cm) lying 
below the plane the liquid nitrogen level. thermo- 
couple bound the outside the finger showed change temperature 
while immersed the liquid. The temperature rose 20° per millimeter 
immediately above the liquid surface, which indicates sharp change the 
temperature the pyrex the liquid surface. Standard adsorbing areas 
cm? and cm? were used the measurements. 


IV. RESULTS 


Figure gives the results six runs with either microamperes 
microamperes, giving gauge pumping speeds differing factor and 


N 

a 

a 

©O 

D 3 

a 

Ww + T 7 


PRESSURE ABOVE ADSORBED LAYER (mmHg 


Fic. isotherm measurements for nitrogen pyrex 77.4° 


with adsorbing areas Most runs were made successive days 
with the adsorbent kept room temperature overnight pressure the 
range. One run (open circles) was performed with pressure 
decreasing, that is, the opening valve (Fig. reduced rather than 
increased the pressure The dashed line Fig. has been drawn visually. 
Monolayer coverage based nitrogen molecular area 16.2 has been 
included for reference. 

The approach physical adsorption equilibrium shown Fig. The 
solid lines show the observed time dependence the pressure following closure 
valve (Fig. three well-separated pressures, for run with increasing 
pressures. pumping gauge were present the system would expected 
that these pressures would become time-independent when physical adsorp- 
tion equilibrium had been achieved. But from eq. (3) continuously 
falls since continuously increases (eq. (6)). Thus pressure falling with 
time expected even when physical adsorption equilibrium has been 
essentially achieved. approximation the expected fall given the 
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OBSERVED PRESSURE (mm Hg) 


TIME IN MINUTES 
AFTER CLOSURE OF VALVE A 


Fic. equilibrium: solid line, observed; dotted line, calculated. 


dashed lines Fig. which have been obtained calculating from 
eqs. (3), (4), and (6) for the runs shown, and obtaining the appropriate pressure 
from Fig. Identity slope between the solid and dashed lines indica- 
tion that physical adsorption equilibrium has been achieved. The dashed line 
has arbitrarily been joined the solid line the end point. The approach 
final equilibrium apparently slow and true equilibrium was probably not 
achieved these measurements. may noted that the true equilibrium 
values are expected lower than those observed. Failure standardize 
the time between closure and desorption (see Fig. contributes the 
scatter Fig. 
DISCUSSION 


Figure compares the present results with some other results the physical 
adsorption nitrogen temperatures near the normal boiling point nitro- 
gen. The units are those commonly used the adsorption literature. 
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Fic. Comparison experimental results the physical adsorption nitrogen near 


Lopez-Gonzalez, Carpenter, and Deitz (1955), Char 27, 77.7° Jura and Harkins 
(1944), 7.77° Bowers (1953), 77.4° Kaganer (1957), aluminum oxide, 
77.4° Kaganer (1958), cocoanut charcoal, 77.6° Law (1955), 77° 
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absolute calibration was made the present experiments and 
the assumption equality between geometric and adsorbing areas was 
arbitrary. Thus detailed quantitative comparisons between the present results 
and those others are not significant but comparisons the general forms 
the results should significant. 

There general continuity between our results for small (approximately 
cm?) plane surfaces and the results Lopez-Gonzalez al. (1955) using 
surface area Law (1955), using small surface germanium 
(approximately finds the first layer complete pressure 
indicate the adsorption quantity corresponding the first layer 
approximately However, despite the absolute difference, the results 
Law (1955) and Kaganer (1957) have the same general form, and, particular, 
Henry’s law proportionality between quantity adsorbed and pressure above 
the adsorbed layer either achieved approached the saturation range 
whereas the present results Henry’s law not valid 
for seeking analytic expression for our results tested 
the equation (Dubinin and Radushkevich 1947) 
employed Kaganer (1957, 1958) for This equation 


the amount adsorbed, 
the pressure above the adsorbed layer, 
the equilibrium vapor pressure above the liquid the same 
temperature, 
dependent the square the absolute temperature, 
and 


where the quantity adsorbed the first layer. Points from the dashed 


line Fig. are plotted Fig. the co-ordinates suggested eq. (7), with 


logio 


p 2 


Fic. Test the equation. 


a 
16 7 
144 
| | 
13 7 
12 | 
ie) 40 80 120 160 
ait — — 


1112 CANADIAN JOURNAL PHYSICS. VOL. 37, 1959 


expressed molecules/cm?. The result straight line for 
This result was not anticipated, particularly view Kaganer’s 
rejection this equation for The value obtained from eq. 
(8) extrapolating the result Fig. Po, done Kaganer 
(1957), molecules/cm?. This almost factor above the 
expected monolayer coverage molecules/cm?, but the discrepancy 
appears within the possible error the experiment. The value 
obtained from the slope Fig. 0.0212, which within the range found 
Kaganer (1958) for various adsorbents (0.0169 0.0564). Thus the application 
eq. (7) our data gives results generally compatible with those Kaganer 
(1957) but much lower saturation values. 

The approach equilibrium shown Fig. while slow, nevertheless 
much faster than the approach equilibrium shown Lopez-Gonzalez al. 
(1955) their Fig. the latter case single isotherm point, taken 
hour after the beginning the process equilibrium, would have yielded 
pressures over the adsorbed layer some three orders magnitude above the 
equilibrium level. The shorter times the present experiment are probably 
the result the thermal contact between bath and adsorbing surface. The 
characteristic thermal time pyrex the experiment 


(9) seconds 


where the thickness the wall, 
the specific heat, 
the density, 
the thermal conductivity. 
For pyrex 77° order few seconds, which experimentally con- 
firmed the small value shown Fig. 


VI. CONCLUSIONS 
Ultrahigh-vacuum techniques are applicable measurements the 
physical adsorption nitrogen small surface areas. 
Use these techniques can reduce the time required reach equili- 
brium. 
The physical adsorption nitrogen pyrex 77.4° was found 
follow the equation for from 
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DESIGN GRADIENTLESS FURNACES! 


ABSTRACT 


Temperature distributions cylindrical furnaces are discussed, and design 
data are given which permit the construction short furnaces which the 
temperature varies less than 0.5% along the central half the length the 
heater. Formulae are also given for the calculation temperature distributions 
conventional furnaces with uniformly wound heaters. 


INTRODUCTION 


Many high-temperature experiments require the absence temperature 
gradients along the working space, and preparing for such experiments 
one faced with the problem designing furnace meet that require- 
ment. Yet information pertinent the design high-temperature furnaces 
exceedingly vague. One sometimes encounters estimates temperature 
distribution particular furnaces, meets such enlightening statements 
“to minimize the temperature gradient inside the furnace, the winding 
density was increased towards the ends the heater’’. However, detailed 
information facilitating the design furnaces meet particular requirements 
seems missing. was alleviate this deficiency information that 
have undertaken extensive study small cylindrical laboratory fur- 
naces. The object view was twofold: primarily, obtain definite design 
formulae for construction furnaces with extended gradientless regions; and 
secondarily, obtain formulae which would permit the calculation tem- 
perature distributions conventional furnaces with uniformly wound heaters. 


DESIGN GRADIENTLESS FURNACES 


2.1 Description the Problem 

The furnace considered here described the following parameters (Fig. 
notice that refers the outside dimensions the tube which the heater, 
length Ly, wound. The furnace filled with insulation conductivity 
K,. water-cooled along and its top and bottom are characterized 
surface conductance h,. The inside tube the furnace contains the 
that term covering all the materials and objects that are inserted into the 
furnace. 

the normal procedure designing furnace for experiment, the 
charge given, and determines both the gradients that can tolerated 
and the dimensions the region over which the gradients have con- 
trolled. Gradientless regions can always obtained furnaces with uniformly 
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Fic. Cross section the cylindrical furnace showing the parameters that describe it. 


wound heaters these furnaces are made very long and thin. Frequently, 
however, such furnaces are unacceptable because the space they occupy, 
awkwardness handling, their large power consumption. Short furnaces 
with extensive gradientless regions can constructed the winding distri- 
bution along the length the heater non-uniform, the non-uniformity 
being function both furnace parameters and charge. One can look this 
non-uniformity being determined the over-all heat flow pattern, which 
can divided conveniently into two parts: the heat flowing predominantly 
through the furnace proper, and the heat flowing predominantly down the 
central tube and charge. These two flow patterns are, course, not indepen- 
dent each other, but turns out that they can treated such good 

The heat flow pattern and temperature distribution the furnace proper 
are functions three geometrical parameters, and the 
winding distribution along the heater, and the ratio heat transfer para- 
meters K,/Lh,. formal analytical solution can obtained for these functions, 
being for the general case double infinite series Bessel functions the 
first and second kind, real and imaginary arguments. This solution much 
too complicated any practical use; was therefore decided tackle 
that part the problem means electrical analogue. 

The heat flowing down the central tube function its geometrical 
parameters and the heat conductivity that tube and the charge. 
was found that this part the heat loss can calculated accurately enough 
the linear flow approximation and subsequently added the main flow 
pattern. The winding distribution the heater then arranged 
meet both heat flow requirements. 
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2.2 Heat Flow through the Furnace Proper 

The flow patterns the furnace proper were determined means 
electrical analogue, the chief characteristics which are here briefly. de- 
scribed. Because the cylindrical symmetry the furnace, only two- 
dimensional resistor network was required. The longitudinal grid the 
network was linear, each subdivision corresponding length L/15 the 
furnace. The radial grid was logarithmic, the junctions being concentrated 
the region interest close had six subdivisions the distance 
The resistors used constructing the network were 0.2% accuracy. 
Power was supplied batteries and the voltages were measured with 
precision vacuum tube voltmeter. The resistors employed the network 
determined the geometrical ratios and and for each set these 
ratios number boundary conditions was applied. These corresponded 
various values Ly/L, K,/h,L, and various power input distributions 
along the heater. The results interest the design gradientless furnaces 
are those which give the power distribution along the heater for constant 
are listed Tables and Table gives the power input factor, from 


TABLE 


Power input factor, for gradientless furnaces 


4.54 3.45 2.24 
4.87 3.47 2.26 
4.56 3.66 2.45 
5.18 3.72 2.47 
4.56 3.91 2.87 
7.07 4.38 2.98 
3.02 2.57 1.85 
4.47 2.84 1.92 
4.56 3.66 2.45 
5.18 3.72 2.47 
9.10 6.85 4.43 
9.53 6.85 4.43 


which the power, required raise the furnace temperature can 
calculated means the equation 


(1) 


where the conductivity the insulation and the length the 
furnace, is, course, dimensionless. Table gives the power distribution 
factor, from which the power distribution along the heater can calcu- 
lated from 


must remembered that the analogue has finite grid, the power 
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TABLE 
Power distribution factor, for gradientless furnaces 


Ly/L = 1 Ly/L = 11/15 La/L = 7/15 


L/ra Kp/hyLao a1 a2 @s as as a6 a7 ao as a4 ao ai ae a3 


4 32 « o @ O08 06 0 0 0 0 0 0 0 3 14 199 0 0O 190 
4 32 0 o 6 G@ 6 4 12 41 468 0 0 4 9 2 20 0 4 5 193 
4 16 oo o & & @ 0 0 0 0 0 Ol 2 86 424 0 5 62 414 
4 16 0 0 0 0 6 14 63 221 624 0 O 11 2 93 474 O 5 62 420 
4 8 © 0 0 0 90 0 0 0 0 O 11 26 59 150 579 O 26 132 610 
4 8 0 6 20 59 141 278 720 2480 0O 13 47 114 280 942 0 27 151 666 
8 16 co) & 8 @ 9g 0 0 0 0 O 8 20 51 142 501 O 26 110 513 
8 16 0 O 8 22 56 133 258 592 2019 O 11 33 93 233 781 O 28 126 554 
+ 16 ro) 0 0 0 0 0 0 0 0 0 Ol 2 86 424 0 5 62 414 
} 16 0 0 0 0 6 14 63 221 629 0 011 2 93 474 O 5 62 420 
2 16 0 0 0 0 0 0 0 & 0 9 145 9 145 
2 16 0 0 0 0 90 0 0 20 358 0 0 0 0 9 151 0O O 9 150 


distribution given not continuous function, but discontinuous values 
pertinent basic interval length For our analogue, and these 
intervals are numbered that refers the middle section the furnace. 
Thus the power dissipated the section the heater counted 
from the center (the central section being and, the center 
plane symmetry, will the same for either end. 

Both the q’s and the a’s are given only for two values the heat transfer 
parameter K,/h,L: zero and infinity. corresponds perfect 
insulation the ends, while perfect water-cooling, which can 
absorb any amount heat without increase temperature. These form 
the limits for any realistic furnace ends, and for most cases the values 
and for K,/h,L equal infinity are close those for equal zero, 
there seemed little point including intermediate values the tables. Some 
practical value. They are included, however, facilitate graphical inter- 

mentioned before, the analogue had grid 6X15. was decided 
run check calculation finer grid, observe the possible effect grid 
size the results. The calculation was done the relaxation procedure 
19/25, and K,/h,L and the results were compared with those inter- 
polated from the analogue tables. The agreement quite remarkable. The 
power input factors agree within 0.1% and the power input distributions, 
shown Fig. also agree very well, indicating that the grid selected for 
the analogue sufficiently fine. 

designing the analogue the conductivity the insulation filling the 
furnace which the analogue represents was assumed independent 
temperature. This is, however, real limitation, the results can easily 
adjusted include insulation whose conductivity not constant. 
using Kirchhoff’s substitution (Carslaw and Jaeger 1959) can very easily 
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Distance along the heater in lL 


Fic. Power distribution along the heater gradientless furnace 16, 


Solid line, given analogue grid broken line, given analogue grid 


shown that for the use insulating powder with con- 
ductivity varying arbitrarily with temperature will leave the power distri- 
bution factor unaffected. The power input will changed only the 
substitution effective conductivity, for eq. where 

1 To 


0 


K,(T) the temperature-dependent conductivity the insulation. must 
stressed that these results are valid only For any 
other K,/h,L the power distribution factor will affected, but, before, 
the values for and form the limits its variation. 


2.3 Central Tube Losses 

The analogue considers only the heat flowing through the furnace proper. 
Heat, course, flows also along the central tube and charge towards the 
ends. This heat flow part the over-all flow pattern, but, mentioned 
before, will treated here independent the furnace body flow, and 
the linear approximation. Such approximation justified because the heat 
flowing down the tube usually only small part the total heat loss. Thus, 
the quantity heat flowing down the central tube, given 


1.80 
| 
140 
1.20 
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where the temperature the center the furnace, 
the temperature the end the furnace, 
the cross-sectional area the central tube, 
the equivalent conductivity tube and charge taken 
temperature 
the distance from the place the heater where provision 
for the supply made the end the furnace. 
the heater arranged that provided its end sections, then 


AX 8/15L—1/2Ly. 


order able use eq. one must make some sort estimate 
our experiments have found very slowly varying function 
being for most furnaces about 180° for 500° and 250° for 
1300° 200° gives close agreement with experimental results. 

Inspection eqs. and reveals bad feature gradientless furnaces, 
namely that the furnace operate over range temperatures, the 
heater winding has made two parts. This arises from the fact that 
the temperature dependence and are quite different, making the 
ratio function temperature. If, therefore, the end losses 
are provided merely increasing the winding density the heater 
its end sections, the furnace will strictly gradientless one temperature 
only. have the furnace operate over range temperatures, separately 
controlled heater must provided supply This auxiliary heater can 
arranged various ways; ourselves have tried three different methods. 
One was insert small end heaters the shape plugs into the central 
tube; another was wind the auxiliary heater top the two end sections 
the main heater; and the third was wind the main heater extending 
only for length and wind the two end sections separately. The 
two end heaters supply then not only but also the power the 
main heater end section. these three arrangements the first gave the best 
temperature distribution, but the third was easiest construct and most 
convenient handle. 


Comparison Design Data with Experiment 

Several furnaces were constructed check the design data obtained from 
the analogue and calculations. The first experiment was made check the 
accuracy the temperature distribution predicted the analogue. This was 
done comparing analogue and observed temperature distributions for 
uniformly wound heater. The results are shown Fig. They agree well, 
the largest deviation being 1.3%. 

Next, gradientless furnace was constructed with parameters 
16, Ly/L 11/15, and length 20cm. The heater was wound 
two parts, with the main winding arranged according Table II, and 
auxiliary winding put over the end sections the main heater. Figure 
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1.00 
1.00 
0.90 
0.90 
0.80 
0.80 
Distance along central tube,in units of Lt Distance along central tube, in units of | 


Fic. Temperature distributions furnace with uniformly wound heater with para- 

Dots, observed distribution; solid line, predicted analogue. 

Fic. Observed temperature distributions for furnaces with parameters 

Solid lines, proper gradientless furnace with two-winding heater; broken lines, uniformly 
wound furnace. 


shows the central temperature distributions observed and 
1000° each sets temperature distributions were measured 
for various ratios The results illustrated Fig. were the best 
obtained and the experimental values and for these distributions 
are compared Table with those calculated from eqs. and Figure 


TABLE III 


Observed and calculated power inputs for gradientless furnaces 


Qn, watts Qr/Qm 
Design temp. 
Calculated Observed Calculated 
600 141 135 208 


1000 277 265 


shows also temperature distributions for furnace with the same parameters 
but uniformly wound heater. 

The winding the heater for the furnace described above was quite compli- 
cated due the variation the winding density along its length. was 
decided test experimentally simplified winding scheme which the main 
heater was uniformly wound along all but the very end sections, where the 
winding density was increased according Table These end sections had 
further winding superimposed them provide The furnace had the 


€ 
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same parameters above. Figure shows the results obtained two central 
the lower temperature the simplification the winding has significant 
effect the temperature distribution. the higher temperature, however, 


1.00 
Wy, 
1.00 
100 
500°C 
0.90 0.90 
(0) ! 2 3 4 5 6 (0) I 2 3 4 5 
Distance along central tube, in units of t Distance along central tube, in units of t 


Fic. Temperature distributions for furnaces with parameters 16, 


11/15. 
Solid lines, simplified gradientless furnace with two-winding heater; broken lines, uni- 


formly wound furnace. 
Fic. Temperature distributions for furnace with parameters 16, 


Solid lines, gradientless furnace with single-winding heater designed for 1250° broken 


lines, uniformly wound furnace. 


impossible obtain true gradientless region with the simplified winding. 
may mentioned here that the values and calculated from 
eqs. and agreed well with those observed for optimum temperature 
distribution. 

The two previously described experiments were carried out with furnaces 
which the heater had double winding. conclude the experiments, 
furnace was constructed with single-winding heater, designed operate 
the temperature 1250° The furnace had parameters 16, 
Ly/L 7/15, and was 20cm long. Figure shows the temperature distri- 
bution obtained for central temperature 500° 1000° and 1250° 
and Table gives the length the region which the temperature varied 
less than 0.5%. seen from that table that the furnace could success- 
fully used the temperature range 1000° 1300° (maximum tempera- 
ture permitted for heater wire), but that the uniform temperature region 


| | 
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TABLE 


Temperatures, power consumption, and length gradientless 
region the furnace with parameters L/ra 16, 


7/15 
Power input, gradientless region, 
watts Calculated Observed 
500 506 
173 1000 1005 
216 1250 1243 


substantially restricted 500° Table also shows the comparison between 
the design temperature, and the central temperature reached with power 
input calculated for that design temperature from eqs. and 


2.5 Discussion and Summary 

practical work the term does not mean precisely zero 
gradient, but usually defined with respect certain maximum limit 
the gradients that can tolerated the experiment. our work 
furnaces, have arbitrarily taken region the maximum 
variation temperature encountered that region equal less than 
0.5%, the temperature being expressed degrees Celsius. With reference 
that definition, all the furnaces that have constructed according the 
above outlined procedures had gradientless region 60% the 
heater length, roughly five times longer than that which would obtained 
with uniformly wound heaters. Although such furnaces are substantial 
improvement the conventional ones, they still fall short the ideal 
which the gradientless region extends over the whole length the heater. This 
shortcoming, believe, mainly due the method employed supplying 
the end losses, Ideally, should supplied wafer-thin circular 
heater located inside the central tube the ends the main heater. 
seems likely that for best results this heater would have wound non- 
uniformly, producing more heat the edges than the center. With such 
means supplying the length the gradientless region should approach 
Ly, its limit. debatable, however, whether the extra extension the 
length the gradientless region worth the trouble constructing such 
complicated end heaters. 

way summary the steps necessary the design gradientless 
furnace will briefly described. The starting point the charge, which 
determined the experiment undertaken. Its dimensions give and 
Ly, the latter the assumption that the gradientless region extends over 
one half the length the heater. Then the ratios and are deter- 
mined, usually compromise between keeping minimum for ease 
handling and maximum keep the power consumption low, and 
the same time avoiding wildly varying winding distribution which com- 
plicates the construction the heater. and are then computed from 
the known parameters, the design temperature and the conductivities 
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the materials employed construction. that connection two references 
may given for the thermal conductivity ceramics; Kingery al. (1954) 
gives extensive data for solid ceramics while Laubitz (1959) discusses the 
conductivity powders. Knowing and with its winding distribution, 
the heater can designed taking into account the nature the power 
supply available and the permissible surface power load the heater wire. 


TEMPERATURE DISTRIBUTION FURNACES WITH 
UNIFORMLY WOUND HEATERS 
3.1 The Pure Radial Heat Flow 

most cylindrical furnaces, the heater wound core whose con- 
ductivity much larger than that the insulation used the body the 
furnace, and the conductivity the charge not negligible, frequently 
happens that the effective conductivity the central tube and charge 
several times that the insulation. Given such case one would expect 
able calculate the central temperature distribution method analogous 
one employed for rod cooling air. fact this type calculation gives 
remarkably good agreement with experimental results for furnaces with 
uniformly wound heaters not more than 10% the heat loss the furnace 
proper occurs through the top and bottom ends. 

Figure shows lines obtained from the electrical analogue correlating sets 
ratios and for which the end losses the furnace proper are 
equal 10% the total loss. There are three lines given, for different values 
Ly/L. These lines apply for the end condition K,/h,L 
gives course end losses, and the end losses for intermediate K,/h,L will 
somewhere between these two limits. All furnaces with uniformly wound 
heaters whose parameters are such give points the left the 


64} 


Fic. Lines correlating sets ratios and for which the end heat losses the 
furnace proper equal 10% the total heat loss. 
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corresponding 10% line have, therefore, end losses less than 10% the total 
loss, and their central temperature distributions can obtained from the 
following equations: 


(4) 
where 


and the distance along the central tube with the origin the center 
the furnace. 
The rest the symbols are previously defined. The power required 
raise the furnace the temperature given 


3.2 Comparison with Experiment 

Two furnaces were constructed check eqs. and one with parameters 
26, 11/15, length 40cm, and the other with 


LOOT 


0.90 
0.90 
0.80 
0.80 


Fic. Temperature distributions for furnace with parameters 16, 
11/15 and uniformly wound heater. 

Dots, observed distribution; lines, calculated distribution. 

Fic. Temperature distributions for furnace with parameters 
11/15 and uniformly wound heater. 

Dots, observed distribution; lines, calculated 
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effective conductivity tube and charge about ten times that the insula- 
tion. Figures and show the experimental and calculated temperature 
distribution while Table compares the observed and calculated power 


TABLE 
Observed and calculated power inputs for furnaces with uniformly wound 
heaters 
Furnace parameters Power input, watts 
11/15 500 111 100 
11/15 1000 278 273 
11/15 500 430 430 
11/15 1000 1318 1321 


inputs required obtain given central temperature The agreement 
between the calculated and observed values certainly remarkable, con- 
sidering the approximations involved eqs. and Both these equations 
were derived not only under the assumption purely radial heat flow, but 
also under the further assumption that and are temperature inde- 
pendent. There simple way which these equations can generalized 
include temperature-dependent conductivities, and the calculation done 
for the comparison shown above average values were used for both and 
Kg, calculated means the formula 


To 0 


should pointed out that the agreement between calculated and 
observed temperature distributions not good towards the ends the 
furnace the central regions illustrated Figs. and the region 
from the observed ones. However, this region the furnace little 
importance, and effort was made pin down this discrepancy. 


SUMMARY 

Design data were presented tabular form which permit the construction 
short cylindrical furnaces with extended gradientless regions. These regions, 
where the temperature varies less than 0.5%, extend over half the length 
the heater and are roughly five times longer than equivalent furnaces 
with uniformly wound heaters. 

Formulae are also given calculate the temperature distribution and 
power requirements conventional furnaces with uniformly wound heaters. 
These formulae are not universally applicable and their limits validity are 
specified. 
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THE DECAY 


ABSTRACT 


has been observed that decays mainly and L-electron capture 
the 2.615-Mev first-excited level Decay known levels above 
the 2.615-Mev level has been shown experimentally not occur more than 
the disintegrations, agreement with analysis the mass 
this region which indicates that such decay modes are energetically 
impossible. Positrons are not emitted more than disintegrations. 
The value the ratio X-rays 2.615-Mev has been measured 
Assuming that the ground state and that the theory 
Brysk and Rose accurate the region near the K-capture threshold, this 
ratio measurement leads value 2.807 Mev for the total decay energy 
This agrees with the value 2.75+0.10 Mev obtained from reaction cycle 
analysis. New measurements the half-life have revealed error the 
previous measurement, the best value now being years within 
factor three. Incidental the main experiment upper limit 0.5% has 
been measured for the positron decay 


INTRODUCTION 
earlier paper (Roy al. 1958) evidence was presented for the pro- 
duction long-lived the fast neutron reaction was 
detected observing the 2.6-Mev y-ray which follows its decay the first- 
excited state Figure shows the level scheme derived 


va 
3.961 Mev 
3.709 
3.475 
3.198 
Mev 
2.703 
72.615 


Ph 


Fic. Level scheme showing known decay and decay observed 
this experiment. y-Ray transitions shown heavy lines are those searched for coinci- 
dence with 2.615-Mev y-ray. Relative intensity (a) from unpublished work Elliott and 


Wolfson. 
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from studies the decay (cf. Strominger al. 1958) and includes 
data both from the work Roy (1958) and the present work. 
The investigation described this paper was undertaken extend the previous 
work, since accurate knowledge the properties could provide 
useful for resolving the mass difference discrepancies that still exist 
the region. 


EXPERIMENTAL 


(a) Coincidence Spectrometer 

From Fig. clear that the decay there should X-rays 
coincidence with the 2.615-Mev y-ray and, sufficient energy 
available, there might also coincident y-rays arising from decay through 
the higher-energy levels 

The coincidence y-ray spectrometer shown diagrammatically Fig. was 
used investigate these possibilities, though some experiments 1}-in. 
diameter }-in. thick crystal replaced the 2-in. 2-in. crystal. This smaller 
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P.H.A. 
E.- A. 
RECORDER 


Fic. Schematic diagram coincidence crystal spectrometer used. 


crystal was both and X-rays and the geometric efficiency 
was computed 49% for thin source placed close the crystal face. 
The detectors were enclosed about inches lead shielding and similar 
thickness lead shielded the small crystal from the stabilizing source. 
The spectrum from the smaller crystal was analyzed the 20-channel 
pulse height analyzer which was gated coincident 2.615-Mev y-ray pulses 
selected from the spectrum the larger crystal. 

The various electronic components were conventional except for the pulse 
height stabilizer which general description given Graham al. 
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(1958, 1086). This stabilizer acted maintain constant pulse height for 
the 662-kev photopeak from the source, and was necessary since, with 
the weak sources available, long counting periods were required and the 
large crystal counter had, itself, inadequate long-term stability. 

Sources and were used calibrate the energy scale 
the smaller crystal spectrometer. Furthermore our bismuth source con- 
tained about ten times much activity produced (n, 
reaction (Eastwood and Roy 1959). Since also decays electron capture 
equilibrium with long-lived precursor was used set the 5-in. 
crystal spectrometer the 2.615-Mev y-ray photopeak. 


(b) Positron Detector 

The standard technique detecting positrons counting coincidences 
annihilation quanta two counters placed 180° the positron absorber 
was used. The apparatus shown schematically Fig. The source was 
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H.T. 
SUPPLY 


SUPPLY 


LINEAR SINGLE SINGLE LINEAR 
CHANNEL CHANNEL 
AMPLIFIER AMPLIFIER 


COINCIDENCE 
MIXER 
SCALER 


Fic. Schematic diagram positron-detecting apparatus used. 


surrounded enough aluminum absorb 2.1-Mev positrons and the whole 
detector setup was enclosed about inches lead shielding. The single- 
channel pulse height analyzers the two photomultiplier outputs were set 
accept only the photopeak from annihilation radiation. The 
efficiency for detection annihilation radiation quanta was about 
whereas the efficiency for random coincidences was The absolute effi- 
ciency for counting positrons was obtained means standardized 
source. The total chance and background coincidence rate was determined 
displacing the source from the axis the two detectors. 


(c) Sources 

The method used prepare radiochemically pure bismuth metallic form 
has been described previously (Roy 1958). For our study the and 
X-ray spectrum with the coincidence spectrometer thin source was 
required order minimize self-absorption effects. The bismuth metal 
powder (120 mg) was therefore spread uniformly over area 7.9 cm? and 
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sandwiched between two films Scotch tape. The superficial thickness the 
bismuth was 15.2 per cm? but because the metal was powder form the 
effective thickness was greater than this. The thickness the Scotch tape 
was per The disintegration rate was obtained from the 
counting rate 2.615-Mev y-rays crystal y-ray spectrometer 
known efficiency (Eastwood and Roy 1959). Assuming that all disintegra- 
tions give y-ray this energy, this rate was disintegrations per second. 
The source also contained about disintegrations per second and 
about y-ray quanta per second (either 270 310 kev) from 
(Golenetskii a/. 1958). 

The source used search for possible positron branching consisted 
100 the same bismuth metal powder enclosed cylindrical polythene 
capsule with walls 1/32 in. thick, which, turn, was placed aluminum 
holder with walls sufficient thickness 1100 per cm?) 
absorb 2.1-Mev positrons. 

RESULTS 
(a) Search for and X-Rays 

Using the coincidence spectrometer described the preceding section 
was found that both and X-rays lead are coincidence with the 
2.615-Mev y-ray Typical spectra are shown Figs. and 

Several checks were carried out confirm that the observed and y-rays 
did fact originate from the decay Despite the fact that the 
and X-rays from the the source were about ten times intense 
those from the the total source activity was low that chance coinci- 
dences these X-rays with the 2.615-Mev y-ray were negligible, amounting 
0.5% and 0.1% for the and X-ray peaks respectively. Chance coinci- 
dences with the similar number higher-energy y-rays from would 
expected least order magnitude less than these. 

separate experiment was carried out show that these coincident X-rays 
arose from the decay and not from lead atoms the shielding which 
had been ionized the K-shell y-ray background particle which 
the same time deposited sufficient energy the large crystal open the 
electronic gate. For this experiment the source was enclosed cadmium 
container in. thick which was sufficient attenuate the X-rays lead 
factor about 1000, and was found that the coincidence X-ray 
counting rate was reduced flat distribution about twice background 
intensity. was therefore concluded that the great majority the coincident 
X-rays originated the source. 

mentioned previously (Roy 1958) y-rays from the known higher 
excited states were found coincidence with the 2.615-Mev y-ray 
(to accuracy This not only ruled out contamination 
the source these high-energy y-rays, but also showed that, the same 
accuracy, decay these higher states takes place. This gives some 
confirmation the estimate from reaction cycle analysis that such modes 
decay are energetically impossible (see Discussion below). 
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RELATIVE NUMBER OF COUNTS 
> 


P.H.A. CHANNEL NUMBER 


TOTAL OBSERVED COUNTS 


diodes 
12 
P.H.A. CHANNEL NUMBER 


Fic. X-rays lead attributed decay observed 2-in. diameter 


(a) X-ray calibration peak obtained from activity source. Curve normalized 


curve (b) peak. 

(b) Observed points and fitted curve for high-energy X-rays counted coincidence with 
2.615-Mev Background has been subtracted and curve normalized equal counting 
time with curve Fig. 

(c) Coincidence background for curve 

(Zero abscissa does not correspond zero pulse height.) 

Fic. Spectrum X-rays lead attributed decay observed diameter 


crystal. Curves show: 

(a) X-ray calibration peak obtained from activity source. Background has 
been subtracted and curve normalized curve peak. 

(b) Observed points and fitted curve for low-energy X-rays counted coincidence with 
2.615-Mev Background has been subtracted. 

(c) Coincidence background for curve 

(Zero abscissa does not correspond zero pulse height.) 


(b) Search for Positron Decay 

positron decay our source could detected. From the counting 
statistics the 99% confidence level was concluded that less than 
positron branching occurred the decay Since the source contained 
about ten times much activity upper limit positron 
branching the decay was set 0.5% this same experiment. 


(c) Measurement X-Rays 2.615-Mev y-Rays 

attempt was made measure directly the electron capture ratio 
order estimate the decay energy from electron capture theory. However, 
self-absorption X-rays the source and lesser extent absorption 
the source-covering and detector window introduced rather uncertain cor- 
rections several hundred per cent these results were without significance. 
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Fortunately measure the ...) electron capture ratio can 
provide essentially the same information about the nucleus does the K/L 
capture ratio. assumed that all decays occur via the 2.615-Mev 
level seems reasonable the light the above observa- 
tions—the ratio X-rays 2.615-Mev will lead directly value 
for the K/(L+M-+ ...) electron capture ratio. Using the crystal spectrometer 
shown Fig. (with crystal) measured the ratio 
X-rays 2.615-Mev y-rays 0.23+0.01. (The necessary corrections for 
counter geometry, self-absorption the source, absorption the crystal 
covering, crystal efficiency, and scattering from the lead-shielding walls have 
been included this value.) 

accurate value 0.957 for the K-fluorescence yield bismuth has been 
determined Mladjenovic and Slatis (1955), who measured the Auger 
electron yield. Combining this with our measured ratio leads value 


(d) Measurement Half-Life 

Absolute intensity measurements have been carried out three different 
samples produced separate pile irradiations These three 
measurements were mutually consistent, but differed about factor 
ten from the previous results (Roy 1958), and check showed that 
calibration error had, indeed, been made the previous work. The cross 
section for the reaction reactor has been re-estimated 
the basis further work other elements (Eastwood and Roy 1959), and 
this new value also differs slightly from that used previously. result 
these two changes now estimated that the half-life 


DISCUSSION 

From these experiments was concluded that decays mainly 
and L-electron capture the first-excited state 2.615 Mev. The 
energy difference between and therefore greater than 2.703 Mev 
(i.e. 2.615 the binding energy Bi). 

value for the exact energy difference may obtained from electron cap- 
ture theory since the K/(L+ ratio indicates that the energy available 
for K-capture not very much larger than the K-electron binding energy 
which case this ratio gives very sensitive measure the decay energy (cf. 
Brysk and Rose 1958, Fig. 6). The theory requires knowledge the type 
electron capture transition which turn may derived from knowledge 
the spins and parities the initial and final nuclear states involved the 
transition. The 2.615-Mev state has been shown experimentally 
state (Elliott 1954). From shell model theory, which well 
confirmed for this particular region (Elliott and Lane 1957), would appear 
that the ground state characteristics will the sum the spins 
the extra proton outside the 82-proton closed shell, and the neutron 
hole the incomplete 126-neutron shell. Thus the ground state will 
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(Nordheim 1951, and also theoretical confirmation de-Shalit 1953) 
most probable that the higher spin state the correct choice. 

assume that the ground state 5+, the electron capture decay 
the first-excited state then first-forbidden unique tran- 
sition (AJ yes) and the theory Brysk and Rose (1958) may used 
obtain relation between the decay energy and the ratio K/L capture 
that completely determinate and independent the form interaction 
the degree parity conservation. the notation Brysk and Rose this 


expression 

order obtain the expression for the K/(L+M-+ ratio for comparison 
with our observations necessary increase the magnitudes terms 
the denominator the above expression factors (defined here 
Which take into account capture from the outer electron shells. 
These factors may obtained good approximation comparing the 
probability density the origin hydrogen-like wave function for the 
and electrons with that for electrons, using Slater screening 
parameters take into account the screening the inner electrons. The 
expression for this factor spin dependent that separate relations are 
obtained for the and the (or terms follows: 


1.272 83) (cf. Wapstra al. 1959, 61), 


(1) K/L 


(3) ra Ly = 3 n° 
(4) (ra = = rn). 


Here the total quantum number and the Slater screening 
parameter (see, for example, Wapstra 1959, 88). Thus for bismuth 
obtain the expression 


which the left-hand side experimentally observed, the and values 
are given Brysk and Rose (1958), and the several values that appear 
may all expressed linear functions any one them making use 
the known electron binding energy differences for the various shells. Solving 
this quartic equation for the neutrino energy A-capture, may 
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then obtain the energy available for the decay the 2.615-Mev 
level and this found 192+3 kev. The error quoted based 
solely the statistics the experiment; however, taking into account the 
plotting and reading errors obtaining the and values from the 
Brysk and Rose curves increases the uncertainty only kev. Adding 
this value the accurately known energy the first-excited state 
obtain the value 2.807 for the total decay energy 

should noted that this value for the decay energy not direct 
measurement but based the assumptions (a) that the ground state 
state, and that the Brysk and Rose theory valid this 
degree precision the energy region where the ratio changing 
rapidly. Since, the authors’ knowledge, there has been experimental 
confirmation that the theory is, indeed, this accurate would unwarranted 
claim have measured the decay energy the precision implied 
the calculations, even assuming that the spin assignment correct. And, 
the other hand, the ground state state, the decay 
first-forbidden transition which case coupling constants and nuclear matrix 
elements appear the relationship corresponding equation (5) above, and 
the known values these quantities are not sufficiently accurate that any 
very useful information may derived from the measurement the 
...) capture ratio. 

interest check whether the spin assignment and the decay 
energy which have derived therefrom lead “‘log value consistent 
with values found for other first-forbidden unique transitions. For unique 
transitions customary include the shape factor (which quite accurately 
known) with the value give value (where the degree 
forbiddenness the (Wu 1955), since found that for these 
transitions the “‘log values form closer grouping that the 
values. For first-forbidden unique electron-capture transitions may 
that 

The lifetime reassessed this paper about years and 
estimated that this accurate within factor three. This gives 
value 11.0+0.5 for which somewhat higher than the range 
values which found for most known first-forbidden unique 
transitions (Wu 1955). the other hand, the values for the three 
other known first-forbidden unique nuclides not exhibit quite 
which our value for has been included for comparison. Most the 
data for this table were obtained from the compilation Robinson and 
Fink (1955); however the log values were computed from equation (5) 


2We are indebted Dr. Lee-Whiting for deriving this relationship well for many 
other valuable discussions and much assistance with the theoretical interpretation our 
results. 
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TABLE 
Properties nuclei that decay first-forbidden unique K-capture 


electron 


Nucleus Decay energy Half-life, years capture 
440+20 100 9.6 
335+10 4.0 2.5 8.4 
(193+4) ~7.5X 105 100 11.0 


*Holland and Lynch (1959). These authors note that this energy value for K“° is inconsistent with that derived 
from the ratio of L-capture to K-capture on the basis of the Brysk and Rose theory, and cite some evidence that 
the decay of K 4° may not be a first-forbidden unique transition. 


above and most cases the uncertainty the value the order 
0.5. From these comparisons conclude that while our “‘log value for 
does not fall within the range normally found for first-forbidden unique 
not inconsistent with the apparent wider spread values 
found for the few cases first-forbidden unique K-capture that have been 
observed. interest note that for the first-forbidden unique 
the 2.615-Mev state the “log value 9.4 which 
also somewhat above the usual range. 

The decay energy may also obtained, principle, from the 
reaction cycle shown Fig. 


CAPTURE 
7-43 
7-40 


Fic. Diagram showing reaction cycle used determining expected decay energy 
(a) Values obtained from charged particle reactions. 
Values obtained assuming 10.65+0.06 Mev. 
(c) Values obtained assuming 10.826+0.020 Mev. 
All energies are million electron volts. Values parentheses are derived from reaction 


cycle analysis (see text). 


Unfortunately some uncertainty exists the neutron-binding energies 
and and value consider two different methods that 
may used obtain these energies. The first method relies charged 
particle experiments which give the limiting values labelled Fig. 
(Harvey 1951, 1953). From these upper limit Mev may 
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derived for the decay energy though the true value may less 
than this limit charged particle group corresponding weak ground 
state interaction has been missed Harvey’s work. This analysis enables 
one set reasonably reliable upper limit 3.00 Mev for the decay energy 
which indicates that the decay via the second higher excited levels 
energetically forbidden (see Fig. 1), was confirmed our 
observations. 

The second method determining the and neutron-binding 
energies makes use photoneutron threshold energies which, unfortunately, 
are usually measured relative the Cu® photoneutron threshold rather than 
absolutely. Thus necessary examine the accuracy the value 
order assess the accuracy other related photoneutron threshold measure- 
ments. 

The most recently derived for the Cu® photoneutron threshold 
energy 10.65+0.06 Mev (Van Patter and Whaling 1957), but find that 
the use neutron-binding energies based this value leads incon- 
sistency the reaction cycle Fig. The binding energies labelled (6) 
Fig. are expressed relative 10.65 Mev for Cu®, the value 7.26+0.05 Mev 
for having been obtained from direct measurements (Van Patter and 
Whaling and the value 4.15+0.10 Mev for having been derived 
from two independent reaction cycle analyses (Hincks and Millar 
These values lead value 2.48+0.11 Mev for the total decay energy 
and this differs from our observed lower limit 2.703 Mev twice 
the standard deviation. 

the other hand, the most recent value for the Cu® neutron-binding 
energy obtained mass spectrometric methods 10.826+0.020 Mev 
(Quisenberry al. 1956). The neutron-binding energy labelled 
Fig. based this value (Chidley al. 1958), and the corresponding 
value for derived previously using this value. These lead 
disintegration energy 2.75+0.10 Mev which agreement with the 
lower limit set our observation X-rays, with the value obtained 
from the Brysk and Rose theory, and with the upper limit set Harvey’s 
work. would therefore appear that the mass spectrometric value for the 
Cu® neutron-binding energy must accepted rather than the value derived 
from betatron and synchrotron experiments are obtain energy 
balance this reaction cycle that consistent with our observations. 


SUMMARY 


summarize, the only observed decay mode and 
electron capture the first-excited state positron emission having 
been detected (to accuracy 5%), nor electron capture higher excited 
states (to accuracy 3%). The half-life has been re- 
measured years within factor three. conclude that 
the decay energy lies between 2.703 and 3.00 Mev, and what appears 


3Note that some adjustments the original values have been made take into account 
the Cu® threshold value used here. 
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the most reliable reaction cycle analysis indicates value 
Mev. This agrees with the value 2.807 derived from theory the basis 
assumed ground state from the value 0.23+0.01 observed for 
the ratio X-rays 2.615-Mev y-rays. Figure shows these conclusions 
diagrammatically. 

by-product this study upper limit 0.5% has been set for the 
occurrence positron decay 
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RADIO STAR SCINTILLATIONS AND 
IONOSPHERIC DISTURBANCES! 


Harr 


ABSTRACT 


The generation mechanism for the ionization irregularities the upper 
atmosphere which are responsible for radio star scintillations considered. The 
general belief that scintillations are related the spread-F phenomenon 
observed ionosonde recordings found inadequate explanation for 
the scintillations Mc/s recorded Ottawa. examination the Ottawa 
recordings shows that there association, both time occurrence 
and geographical location, with those ionospheric disturbances that are usually 
considered due incoming solar particles. Since other workers more 
southerly geomagnetic latitudes have associated their scintillation observations 
with the spread-F phenomenon which they consider independent 
auroral activity, would appear that two mechanisms, least, are responsible 
for the radio star fluctuations: namely, the precipitation solar corpuscles and 
mechanism linked with the spread-F phenomenon. The former seems pre- 
dominate high latitudes, the latter probably predominant low latitudes, 
while both mechanisms probably are operative middle latitudes. 


INTRODUCTION 

Radio star scintillations have been studied for some years, but, while they 
are known result from irregularities ionization the upper atmosphere, 
little agreement has been reached the origin the irregularities. Some 
authors favor explanation involving the region and the phenomenon 
known spread (Dagg 1957a), others would place the irregularities 
the region and link them the sporadic-E phenomenon (Bolton, Slee, and 
Stanley 1953). neither case, however, the generation mechanism for the 
irregular structures understood fully. 

Several recent papers the subject seek locate the mechanism for the 
production the ionospheric irregularities within the terrestrial atmosphere. 
Booker (1958) argues favor atmospheric turbulence height between 
the and regions. However, fails account for the observed nighttime 
maximum scintillation daytime maximum should fit the turbulence 
model much better. Dagg (1957c, reviews number existing theories, 
pointing out their shortcomings, and develops suggestion that scintillation- 
producing F-region irregularities are transferred there from the dynamo region 
where they arise through turbulent winds. Neither these hypotheses seems 
entirely satisfactory since adequate explanation presented for the known 
dependence scintillation ionospheric disturbances high latitudes. 

Dagg actually raises interesting point when considers the F-region 
irregularities that have been observed Peterson (1955). While such 
irregularities may well explain the high-latitude scintillations, Dagg feels they 
cannot explain scintillations observed equatorial regions, and discards 
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the idea. This action seems too drastic since one has yet satisfactorily 
examined the possibility that more than one mechanism responsible for 
scintillations. the latter were the case, one might not expect much similarity 
between the results for high and low latitudes. Indeed, the hypothesis that 
two separate processes are involved would seem account for most the 
observed behavior scintillations and, should one the processes involve 
the phenomenon known spread would also account for the apparent 
close connection with irregular ionization structures the region, under 
certain conditions. Koster (1958) has already pointed out the differences 
which exist between scintillations observed equatorial station and those 
observed higher latitudes. there link between F-region irregularities 
and scintillations one would expect that spread-F conditions low and high 
latitudes would dissimilar, and such has been shown the case since 
spread correlates positively with magnetic activity high latitudes, and 
negatively low latitudes (Lyon, Skinner, and Wright 1958). Consequently 
would seem worth while examine the possibility that scintillations 
high latitudes are primarily due incoming particles which produce disturbed 
ionosphere conditions, while those equatorial regions are primarily due 
some other mechanism, possibly one that contained entirely within the 
ionosphere. Under such circumstances would seem highly likely that 
scintillations temperate latitudes would due combination the 
two effects. (An inherent requirement this thesis that low-latitude 
spread not produced the precipitating solar particles, which con- 
sistent with the observations Lyon, Skinner, and Wright (1958).) 

The phenomenon spread which observed ionosonde recordings 
implies irregularities the F-region ionization. Since this condition 
common one, not all unreasonable expect that radio waves that 
arrive from outside the ionosphere would scattered such F-region 
irregularities. Briggs (1958a) has recently been able show detailed corre- 
lation between spread-F conditions and scintillations for the ionosphere 
the same geographical region, which should remove some doubts the 
matter. However, the lack perfect correlation his data leads one 
suspect that has not yet provided the complete answer, and perhaps 
combination circumstances would have resulted better explanation. 
Briggs suggested that the scintillation data was considering were not 
the type which conceivably could linked with auroral activity, but his 
proviso seems based nothing more than that positive correlation with 
spread-F echoes was found and that the scintillations were 
phenomenon. This does not seem, itself, adequate basis for the selection 
his data since aurora known occur latitudes where the Cambridge 
scintillation data were obtained (Vestine 1944). Vestine’s map shows that 
the line sight the radio star passes through the ionosphere 300-km 
height geographical region for which there occurrence aurora 
upper transit the source, Cassiopeia and about 10% occurrence 
lower transit. Therefore might expected that the contamination the 
scintillation data auroral activity should minor factor. However, this 
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argument weakened the well-known fact that there not one-to-one 
correspondence between visible auroral activity and the ionospheric ionization, 
associated with auroral disturbances, that gives rise the radio reflections 
the very-high-frequency and ultra-high-frequency bands. Collins and For- 
syth (1959) have found considerable scattering very-high-frequency waves 
‘radio-aurora’ subauroral latitudes, even times when visible aurora 
was not present. Accordingly, the Cambridge scintillation data may not 
entirely free that type auroral activity which associated with incoming 
particles. 

should noted that high-latitude spread occurs when the whole 
ionosphere disturbed incoming particles (cf. Hartz and the 
possibility cannot overlooked that some other ionospheric region respon- 
sible for the high-latitude scintillations the same time that the signals 
are spread: when the radio star signal passes through active aurora 
100-km height and throughspread are the scintillations produced 
the region? Certainly, such times the region known possess 
irregular ionization structures that have the shape and motions required 
produce the observed scintillations (cf. Lyon and Kavadas 1958), 
nounced changes the scintillation indices have been observed during 
auroras (Little and Maxwell 1952). 

The Vestine map shows that the auroral activity Ottawa higher than 
Cambridge. The line sight from Ottawa the Cassiopeia source 
passes through the ionosphere 300-km height region which there 
10% occurrence visible aurora upper transit and 50% occurrence 
lower transit. Consequently, the effects auroral activity should more 
noticeable the Ottawa scintillation data. These data will considered 
the following sections, and, particular, dependence particles solar 
origin impinging the ionosphere will sought the data. Since the 
term ‘auroral activity’ rather indefinite, several measures particle preci- 
pitation will considered and correlated with radio star scintillations. This 
comparison indicates that mechanism linked with general ionospheric 
disturbances responsible for the observed fluctuations Ottawa. 


SCINTILLATION MEASUREMENTS 

Radio star scintillations have been recorded Ottawa frequency 
Mc/s for number years. phase-switching interferometer receiver 
more less conventional design (Ryle 1952) has been used monitor the 
intense source Cassiopeia The recording station was located (latitude 
longitude 75.9° W.) that the source never set. Continuous records 
were obtained with time constant 1.5 seconds the meter circuit. The 
charts were analyzed assigning index the amplitude and the rate 
fluctuation for each hour recording. Amplitude indices ranging from 
were chosen the basis the ratio the fluctuation amplitude the 
amplitude the sinusoidal recorder pattern. Examples the respective 
amplitude indices are shown Fig. The rate indices used ranged from 
where meant scintillations, was assigned describe the slowest 
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Samples the scintillation amplitude indices. 
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fluctuations observed, and higher indices for increased scintillation rates with 
characterizing the most rapid rates recorded. Examples the rate indices 
are shown paper Reid (1957). The reader may note that the present 
amplitude indices differ somewhat from those described earlier paper 
(Hartz these new indices permit finer gradation fluctuation ampli- 
tudes and not affect any the basic conclusions the previous analysis. 

The earlier work showed that the scintillations observed Ottawa did not 
agree with results obtained elsewhere (Hartz 1955a). particular, scintilla- 
tions were observed all hours the day, and, while nighttime maximum 
occurrence was noted, was not large maximum. Moreover, definite 
correlation with the spread-F phenomenon could established. Further con- 
siderations disclosed that there was significant angle-of-elevation dependence 
the data, such that low angles elevation the fluctuations occurred 
more frequently. easy see from the geometry involved that the angle 
elevation the source changes the receiving site, the distance from the 
zone maximum auroral occurrence that ionospheric region given 
height through which the rays pass also changes. Consequently, the observed 
angle-of-elevation dependence for scintillation occurrence can equally well 
described dependence the distance that the active region the iono- 
sphere from the so-called Auroral Zone. fact, the latter explanation 
much more acceptable one since Ryan and Harrower (private communication), 
the basis scintillation measurements made Kingston, Ontario, have 
found that the angle-of-elevation dependence for sources the southern sky 
does not correspond with that found the northern sky. 

the Auroral Zone considered that region which most the 
incoming solar particles precipitate, seems worth while explore the con- 
nection between the Ottawa scintillations and the incoming particles. For 
this necessary quite definite the criteria which establish the 
occasions which particles are said entering the ionosphere. Direct 
observations the particles are virtually impossible, that indirect measure- 
ments must used continuous data are obtained. the following 
sections several such indirect measures particle precipitation will 
considered and their relationship radio star scintillations will examined. 


SCINTILLATION RATES 


The data for the rates scintillation for the 2-year period from July, 1955, 
July, 1957, were summed and plotted function solar time, shown 
Fig. this diagram the percentage the time that each fluctuation 
index occurs shown for each hour the day. can seen that the high 
scintillation rates predominate night and the lower rates the daytime. 
This shown better Fig. where the percentage occurrence high 
scintillation rates, that indices and plotted curve For interest, 
the percentage occurrence high scintillation amplitudes (Hartz 
also presented here curve evident that some correspondence exists 
between the two measures both curves indicate higher activity night 
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Fic. occurrence (A) high scintillation rates, and (B) high scintillation 
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than during the daytime. hour-to-hour correlation between the unsmoothed 
amplitude and rate data for the 2-year period yielded correlation coefficient 
+0.44. 

The scintillation rate data were also summed and plotted function 
sidereal time. This presented here Fig. the ordinate shows the per- 


PER CENT 


° 4 8 12 16 20 24 
SIDEREAL TIME (HR) L.T. 


Fic. occurrence the various scintillation rates function sidereal 
time. Lower transit the source taken hours. 


centage occurrence the various scintillation rates and the abscissa the 
time sidereal hours after lower transit the source. may seen that 
low angles elevation the source the higher rates (particularly index 
predominate, while high angles elevation the lower rates predominate. 
(For Ottawa the angles elevation lower and upper transits Cassiopeia 
are respectively 13.9° and 76.9°.) again one considers only the high 
scintillation rates, namely, indices and the graph shown curve 
Fig. obtained, while curve that figure gives the percentage occurrence 
the high scintillation amplitudes (Hartz 1955a) for comparison. Again 
apparent that there considerable correspondence between the high 
scintillation rates and the high scintillation amplitudes. 

From the foregoing readily seen that, the average, there com- 
ponent the scintillation data which varies with the sidereal period 
hours minutes seconds. Whether this can considered due the 
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SIDEREAL TIME (HR) 


Fic. occurrence (A) high scintillation rates, and (B) high scintillation 
amplitude function sidereal time. 


changing path length for the radiation through ionospheric layer, the 
changing distance from the Auroral Zone the ionospheric region where the 
radiation penetrates, not immediately clear from only this evidence. Likewise, 
the scintillation data contain component which varies, general, with 
period hours. doubt this related some way ionospheric con- 
ditions which depend solar radiations. Before any conclusions could 
drawn the relation the radio star scintillations incoming solar particles, 
these respective components had removed from the scintillation data. 
This was done applying 24-hour running average the data. The 
smoothed scintillation rate indices were then examined for evidences 
incoming corpuscles 3-hour basis. The results are presented the 
following sections. 


MEASURES DISTURBANCES RESULTING FROM 
INCOMING PARTICLES 


disturbances may classified according their generation 
mechanisms into two main types: those which are produced abnormal 
electromagnetic radiations, and those which are attributable charged 
particles which enter the ionosphere. This latter category has been shown 
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contain those events which are characterized magnetometer disturbances, 
visual auroral displays, and abnormal variations the electron density 
the upper atmosphere. generally accepted that the responsible particles 
have solar origin, and that the particles enter the earth’s environment 
they produce abnormal conditions the ionosphere widespread and 
complex that they cannot described adequately any way other than 
general ‘ionospheric disturbance’. The accompanying magnetic storm, 
auroral features, etc. are only particular features the general disturbance 
resulting from the precipitating particles. For the present purposes such 
features will used indicators the incoming solar particles. 


(a) Geomagnetic Data 

The planetary indices magnetic activity were used indications 
general ionospheric disturbances that result from incoming particles. was 
felt that the planetary indices provided measure the currents that flow 
the ionosphere, probably even the ring current such exists, conse- 
quence the solar corpuscles. has already been established that there 
positive correlation between scintillation indices and the planetary magnetic 
indices and also between scintillation indices and local magnetic indices, with 
the coefficient the latter being greater than that the former (Costain 
1955; Dagg 1957b). This expected both scintillations and local 
magnetic variations are measures related ionospheric conditions the 
vicinity the observing station. However, for the present analysis the 
relationship the scintillation indices incoming particles was sought and 
the planetary indices were considered more reliable for this than are local 
indices magnetic activity. 

The indices were smoothed applying 24-hour running average 
make them compatible with the scintillation rate data, and then correlation 
was carried out 3-hour basis. For the 2-year period, July 1957, 
correlation coefficient +0.54 was obtained for these two sets data. 


(b) Data 

was difficult decide single ionospheric parameter that could 
easily measured which would give reliable indication incoming particles. 
While certain observations the region means ionosonde are 
probably acceptable most the time, difficulty arises just times 
disturbance, because D-region absorption. choice was finally made 
the transmission quality long-range high-frequency communication path 
which skirted the Auroral Zone. The data were obtained from the Canadian 
Overseas Telecommunications Corporation for the trans-Atlantic circuit 
between Montreal and London. Hourly indices ranging from had been 
assigned according reception conditions were excellent, good, fair, poor, 
very poor, ‘blacked-out’. While, strictly speaking, this was not quanti- 
tative measure, could still serve reliable indication abnormal con- 
ditions since the operator made comparison what would normal 
reception for him. Further, the operator would adjust for the changing 
transmission conditions different times the day and year accordance 
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with the changing ionosphere, and list transmission quality relative what 
would normal for that particular time. 

Because the reflection process was involved, this particular measure provided 
disturbance data for only the lower regions the ionosphere. sure, 
similar data could have been obtained from vertical soundings, but would 
have been less convenient form since the automatic recorder cannot 
make comparison with normal ionosphere conditions can the manual 
operator. Moreover, the trans-Atlantic path was preferred any other 
operational circuit since skirted the Auroral Zone, which considered the 
most likely region for the precipitation incoming charged particles. 

The trans-Atlantic quality data for the 2-year period, July 1957, 
were smoothed applying 24-hour running average. this way, minor 
fluctuations were smoothed out, were the daytime effects the sun’s 
ultraviolet radiation. The major variations that remained the data could 
attributed corpuscular radiation. These smoothed data were then 
correlated with the smoothed data and the smoothed scintillation rate 
indices 3-hour basis. Correlation coefficients +0.66 and +0.29 re- 
spectively were obtained. 

(c) Data 

observer near the Northern Auroral Zone sees the aurora borealis almost 
every clear night, but further south the frequency overhead aurora decreases, 
until southern Canada fairly rare phenomenon. While aurorae are 
considered occur consequence incoming charged particles, quanti- 
tative measure the particles difficult establish from the visible aurora 
since many parameters are involved. possible measure the particles 
seems contained the latitude which the aurora seen, since sightings 
southern latitudes are known occur mainly times intense dis- 
turbance. 

Data the position the aurora function time for recording site 
Saskatoon (latitude 52.1° N., longitude 106.6° W.) were made available 
the author Dr. Hunten. From these, the most southerly extent 
the aurora for each night observation was chosen measure the 
auroral disturbance. The lower frame Fig. shows these observations 
terms geographic latitude for the month January, 1956. (Note that the 
latitude range was all that could reported reliably from the one 
site, and occasions when the aurora extended beyond the limits observa- 
bility the south, only the limiting latitude was given.) the diagram, the 
closed circles represent weak display, while the open circle represents 
more intense display. 

The auroral data were discontinuous make the calculation 
correlation coefficient very difficult and none was attempted. However, when 
portrayed graphical form, which Fig. example, was found that 
the maxima the scintillation rate curves tended occur those days 
when aurora was observed the lower latitudes, and was increased 
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TRANS-ATLANTIC QUALITY 
SCINTILLATION RATE 


INDICES 


LATITUDE AURORA 


JANUARY 


ionospheric disturbance for the month January, Plotted here 
are the smoothed trans-Atlantic quality indices ranging from for excellent for ‘blacked- 
out’, the smoothed scintillation rate indices, the smoothed planetary magnetic indices, and 
the southermost extent the aurora borealis each night for which records were available 
(the open circle represents very intense 


Also plotted Fig. for comparison, are the smoothed data trans- 
Atlantic propagation quality, radio star scintillation rates, and planetary 
magnetic indices for January, 1956. may seen that the four different 
measures presented the diagram are not independent, substantial 
degree correlation exists between the respective sets data. particular, 
the large significant disturbances seem portrayed reliably any one 
the four, and only with the lesser disturbances that differences the 
observed effects occur. Such differences can probably attributed differ- 
ences the latitude height the ionosphere which the particular 
measurements apply. the whole, believed that each phenomenon 
plotted here manifestation corpuscular radiation from the sun. 
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The reader reminded that the planetary magnetic indices are indicative 
the total motion charges the ionosphere which arise from the incoming 
corpuscles. abnormalities the electron densities occur the lower iono- 
sphere these would probably affect the magnetic indices well the trans- 
Atlantic quality figures, and hence the correlation coefficient +0.66. The 
motions small regions ionization high levels the ionosphere, which 
probably contribute the observed scintillation indices, also add the 
current fluctuations which determine the values, and hence the correlation 
+0.54. That higher ionospheric regions than those which deter- 
mine the high-frequency quality data are involved some the scintillation 
measurements inferred from the low correlation coefficient between scintilla- 
tion rates and the short-wave quality figures, namely +0.29. this the 
case, might expected that the trans-Atlantic quality data and the 
scintillation rate data were combined, and then compared with the indices, 
higher correlation coefficient would result. When the indices for the 2-year 
period were added together, and the correlation was carried out 3-hour 
basis, the result was correlation coefficient +0.72 between the trans- 
Atlantic quality and scintillation rate combination and the data. This 
result consistent with the general belief that scintillations arise above 
the region the ionosphere. also consistent with the hypothesis that 
part the scintillations arise the region above the level where high- 
frequency absorption occurs, since D-region absorption the major factor 
that affects the trans-Atlantic quality data. 


DISCUSSION 


There can doubt that the radio star scintillations observed Ottawa, 
and probably other stations comparable latitude, are definitely associated 
with the large ionospheric disturbances which result from incoming solar 
particles. (In subsequent paper the author will show the dependence the 
significant ionospheric disturbances this period specific solar events, 
and thereby establish the solar origin the particles each case.) There can 
also doubt that there definite diurnal dependence the scintillation, 
but this should cause surprise since there are other geophysical phenomena 
associated with incoming particles which also exhibit diurnal variation. 
These include local magnetic disturbances, radar reflections from auroral 
ionization, and the scattering very-high-frequency radiations the forward 
direction auroral ionization. fact, the similarity the curves Fig. 
the diurnal variations each these three phenomena sufficiently great 
indicate itself definite relation between the respective sets data 
(cf. Bullough and Kaiser 1955; Collins and Forsyth 1959). Kaiser (1956) has 
already drawn attention this similarity, particularly connection with 
the minimum midnight, and has drawn the obvious conclusion that some 
portion the scintillations least must produced the region the 
intense irregularities which are shown exist the auroral radar echoes. 

There seems little argument against Briggs’s conclusion that the 
scintillations not associated with auroral activity are closely related the 
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F-region irregularities that produce the spread-F echoes ionosonde records 
(Briggs 1958u, But what the scintillations that are associated with 
auroral activity? Here, the matter the height the scintillating region 
the ionosphere requires more intensive consideration. F-region irregularities 
ionization capable scattering very-high-frequency radio waves are known 
exist times when aurora present. The author has observed forward- 
scatter signal Mc/s 3000-km circuit time when aurora occurred 
along the path: for path this length the ionosphere below 170 the 
mid-point was not illuminated the transmitter nor observed the receiver. 
Peterson (1955) and Nakata (1958) have also observed return echoes 
from F-region irregularities which were associated with auroral phenomena. 

the other hand, irregularities capable scattering very-high- 
frequency radiations are well known. Recent work Collins and Forsyth 
(1959) signals scattered the forward direction has shown that several 
different types irregularities are involved, not all which are associated 
with visible auroral phenomena. They were able determine the height the 
scattering region and most cases the results showed heights less than 
100 km. 


the Ottawa area showing lines equal percentage auroral occurrence, after 
Vestine Also shown the locus intersection with the 300-km level the ionosphere 
line joining Ottawa Cassiopeia The numbers the locus represent sidereal time 
hours after lower transit. 
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Apparently, during disturbances that result from incoming solar particles 
all ionospheric heights are affected some degree other. From the informa- 
tion available does not seem likely that any one height can designated 
the exclusive origin the scintillations that occur such times. This 
conclusion can substantiated further through the following argument. 

the ionospheric region responsible for the scintillations could confined 
limited height range, should possible determine this height from 
the observed variation scintillation occurrence with sidereal time. Vestine’s 
map auroral occurrence can considered representative the geograph- 
ical distribution ionospheric structures associated with incoming particles, 
particularly since similar the map which portrays the magnetic dis- 
turbances (Vestine 1944). Fig. shown portion this map covering the 
Ottawa area, and including the locus points intersection with the 300-km 
level the ionosphere lines drawn from Ottawa the source Cassiopeia. 
The figures the locus are the time measured sidereal hours after lower 
transit the source. the percentage occurrence aurora each hour 
obtained interpolation process from this diagram, the values plotted 
Fig. are obtained. will seen that there considerable similarity between 
this curve and the curves Fig. One might expect that somewhat different 
choice height for the scattering layer could made which would permit 
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better agreement with the experimental curves Fig. but this has not been 
possible. Almost any height less than about 600 would have resulted 
curve similar shape, but close fit could only obtained through normali- 
zation procedure which would difficult justify. This may due the 
fact that the treatment has been too elementary and the scattering properties 
the irregularities have not been considered properly, the possibility 
that the irregularities are actually distributed through extensive range 
heights. the light the foregoing argument and F-region structures 
the latter suggestion should given serious consideration. 


SUMMARY 


The radio star scintillations observed Ottawa have been examined for 
dependence incoming charged particles which produce ionospheric dis- 
turbances. The planetary indices magnetic activity, the quality reception 
high-frequency path that skirted the Auroral Zone, and the lowest 
latitude which aurora was observed were taken reliable measure 
ionospheric disturbances caused solar particles. The scintillation data, from 
which diurnal variations had been smoothed, correlated significantly with each 
the three disturbance measures, and therefore must taken repre- 
sentative ionization structures produced the precipitation corpuscles. 
has not been possible determine the height the scintillating region 
the ionosphere directly, but known that during ionospheric disturbances 
man-made transmissions have been scattered the forward direction 
irregularities the region and the region the ionosphere, and therefore 
difficult exclude the possibility that the scintillations are produced 
extensive range heights. 

latitudes where the auroral activity much less the scintillations are 
known correlate with spread-F echoes and presumably originate the 
region. should not taken that such scintillations are not found the 
Ottawa records, but rather that they not predominate. The best explanation 
that the recordings show the presence both types scintillations. 
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THE PHOTODISINTEGRATION NEON! 


ABSTRACT 


proportional counter, filled 4-atmosphere pressure with pure neon, and 
cylindrical gridded ionization chamber, filled 9-atmosphere pressure with 
neon containing some helium, have been irradiated with the gamma-ray flux 
from thick lithium target bombarded with 500-kev protons. Pulse height 
analysis led the assignment the following photodisintegration cross sections 


millibarns: 
0.05 the ground state 0.67 the 6.06- and 6.14-Mev 


levels 1.80 the 6.91- and 7.12-Mev levels 
0.086 the ground state, about 1.0 the 6.06- and 6.14-Mev 


levels; 

the ground and low-lying levels 110, 197 kev 
about the 1.35- and/or 1.57-Mev levels 

0.76 the ground level. 

The cross sections are good accord with the selection rules 
put forward Gell-Mann and Telegdi for the photodisintegration 


even nucleus. 
INTRODUCTION 


study the disintegration neon nuclei mass energetic photons 
provides opportunity examining two aspects present ideas about 
photodisintegration processes. the first place since alpha-particle 
model, may regarded bipyramid alpha-particle groupings, then for 
photon energies not big enough disrupt the alpha structures the charge 
center will not oscillate with respect the mass center under the influence 
electric dipole radiation. Thus the absorption process would expected 
arise from magnetic dipole electric quadrupole interaction. Following 
the excitation the compound nucleus, the binding between the alpha- 
particle configurations relatively small, then the ejection alpha par- 
ticle should quite probable process. The reaction was 
fact the first observed photodisintegration process neon. Erdman and 
Barnes (1953), and greater detail Erdman (1953), studied this reaction 
using the radiation from lithium bombarded protons and observing the 
disintegrations produced neon contained parallel plate gridded ioniza- 
tion chamber. They found surprisingly small cross section for the dis- 
integration proceeding the ground state and noticed that considerably 
more transitions occurred the and 7-Mev excited levels Preston 
(1950) has given theoretical analysis the process which 
might expected apply the photo-alpha disintegration neon yielding 
cross section similar order magnitude. Gell-Mann and Telegdi (1953) 
investigated the consequences the charge independence nuclear forces 
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the absorption photons light nuclei and showed that isotopic spin 
conservation should lead strong, clear-cut selection rules for the case 
nuclei such Again, observation the alpha energies 
the disintegrations via excited states are quite apparent and this 
technique provides another opportunity checking the non-existence, accord- 
ing the alpha-particle model, low-lying states the nucleus. 

second process interest the photoproton disintegration 
Studies the gamma-ray yield from the reactions 
(Willard 1952; Sinclair 1954) show that the levels the com- 
pound nucleus are reasonably well separated even excitations 
Mev. Consequently, the principle reciprocity, the cross section 
could show quite sharply resonant behavior which might detectable 
varying the gamma-ray energy over small range. direct numerical 
check too the applicability the reciprocity relation, 
might also feasible, since with modern van Graaff accelera- 
tors the same excitation can reached either proton bombardment 
fluorine absorption 17.6-Mev quanta The prime difficulty 
with the measurement would, course, lie measuring only the weak 
direct-capture radiation component the large flux other gamma radiations. 

Because neon excellent counter gas, except for its low breakdown 
voltage, the photodisintegration process leading charged particle emission 
may conveniently studied means proportional counter ionization 
chamber technique, pulse amplitude analysis giving the energy released 
the counter from the disintegrations. Moreover, the range alpha particles 
very much less than that protons comparable energy, that change 
the gas pressure permits clear distinction between groups the pulse 
height distribution arising from reactions and those from re- 
actions. This was the method used our study the photodisintegration 


EXPERIMENTAL APPARATUS 

(a) Source Radiation and Monitoring 

The gamma radiation was generated bombarding thick evaporated 
layer metallic lithium with resolved 500-kev energy proton beam 
400 microamperes from the A.N.U. cascade generator. 
This radiation comprises two components, one sharp, energy 17.6 Mev 
with spread only kev, and the other broad, centered around energy 
14.8 Mev spread about Mev. The intensity ratio, 
the beam amounts 1.6/1 (Stearns and McDaniel 1951). The yield, 
according Fowler and Lauritsen (1949), for clean fresh target the 
440-kev resonance about 1.2 gammas per microcoulomb, 
per proton, and this increased about 10% 500 kev the 
contribution from protons above the resonance energy. 

The radiation was monitored with brass Geiger counter similar con- 
struction that calibrated Barnes al. (1952) having efficiency 
the mean effective photon energy 16.4 Mev 2.5 counts per photon per 
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square centimeter and this figure believed reliable within +10%. 
Generally the flux for given charge was only about one-half that computed 
from the yield data; only when the target had been freshly laid down was 
the full yield obtained. 

Background counts with the beam stopping clean copper target showed 
that very few counts the monitor Geiger arose from rays from the 
accelerator, the counting rate remaining about counts/minute with the 
beam off. Carbon build-up the target and consequent radiation 
from the y)N™ reaction, with resonant yield 7.3 gammas 
per proton 456-kev bombarding energy, was largely avoided placing 
liquid air trap immediately front the target and running the target 
quite hot. There was evidence after any run significant 10-minute 
positron activity from build-up. Checks made with neutron-sensitive, 
gamma-insensitive scintillation counter showed that the number neutrons 
originating from the bombarded target assembly was negligible. 


(b) The Proportional Counter 

With the amount very neon available was decided first 
investigate the (y, reaction using proportional counter, conventional 
cylindrical design with sensitive region in. I.D. in. long; the .005-in. 
diameter central wire was provided with shield and field tubes, the latter 
extending inches inside the walls each end. After purifying the gas 
circulation over hot calcium, the counter was operated low gas gain and, 
with differentiating and integrating time constants microseconds, the 
width the 5.16-Mev alpha group from plutonium energy calibration 
source was 6%. 


(c) The Gridded Ionization Chamber 

order study the process was necessary use much 
greater stopping power the counter, i.e. higher pressure and bigger 
volume. Consequently cylindrical geometry, gridded ionization chamber 
was built shown Fig. The main advantage such chamber over 
proportional counter for this work that the field the walls may made 
much larger, bearing mind that the total applied voltage limited 
the breakdown voltage neon, thereby reducing the effects possible con- 
tamination. The cylindrical design makes better use the available gas than 
does parallel plate chamber. 

The neon available fill this chamber was rather impure, containing 
particular nitrogen and helium. This gas was passed about 
through system consisting activated charcoal trap 
liquid air, followed hot calcium purifier, and followed 
further charcoal trap. Mass spectrograph analysis, kindly undertaken Dr. 
Richards this laboratory, showed the final gas composition 
70.05%, Ne” 7.15%, 22.80%, the nitrogen having been successfully 
removed. The presence the helium appeared raise the breakdown potential 
the gas and was possible run with 4000 volts potential difference 
between anode and wall total gas pressure 132 this limit actually 


t 

H 

13 

i 

H 

? 

t 


STHINI 


OOM 1331S 
OVD 


BOONY 


wae Lee 


CANADIAN JOURNAL PHYSICS. 37, 1959 


@ 


wninoinid 


1156 
| 
| 
‘Al 
i 
| 
| 
| 
Eat 


HAY AND WARREN: PHOTODISINTEGRATION NEON 


GaovdS AWWNOA 1S O€ 


1157 
| 
| 


1158 CANADIAN JOURNAL PHYSICS. VOL. 37, 1959 


being set the power supply noise the collector was run high potential. 
The grid potential was set one-half the anode potential indicated 
calculation, and tests made with various grid voltages showed that this choice 
was about optimum, though neither resolution nor pulse height was very 
sensitive function grid voltage. Calculation (Bunemann al. 1949; Wilkin- 
son 1950; Whitehouse and Galbraith 1950) showed that the grid-shielding 
efficiency under these conditions was 0.96 while its transparency electrons 
was unity. The minimum value E/p the chamber was 0.019 volts/cm 
per compared with 0.005 the proportional counter, and with this chamber 
there was much less difficulty with contamination. The calculated maximum 
electron-collection time for 10-Mev alphas range this pressure 
was about microseconds, while for 4.7-Mev protons range about 
was microseconds, taking the electron drift velocity 
per microsecond the gas (English and Hanna 1953); experimental results 
showed that these figures were approximately correct. 

was found possible achieve width for the plutonium alpha group* 
less than for values the amplifier integration and differentiation 
time constants the range microseconds. The wall thickness was 
reduced inch around the sensitive region reduce electron pile up. 
Nevertheless, with the gamma fluxes used this calibration peak was con- 
siderably broadened, about 10% width, the electron noise, and colli- 
mation the beam made little improvement. The group widths obtained 
photodisintegration runs were rather broader still (~12%), indicating small 
additional spread caused either the events being distributed uniformly 
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height spectrum from neon-filled ionization chamber irradiated with 
radiation and background. 


*The plutonium source was kindly supplied the Atomic Energy Research Establishment, 
Harwell, England. 
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throughout the sensitive volume the nature the events themselves. 
This width was obtained with integration time constant microseconds 
and differentiation time microseconds, the latter being short 
was compatible with the variations rise times. The time required for 
experimental run was found set electron pile up, not the total 
gamma flux available. Consequently, some improvement resolution might 
achieved the provision thin low windows for the radiation which 
would also permit the use higher proton beams, i.e. greater than 200 micro- 
amperes, reduce the importance counter background. The natural back- 
ground this chamber was quite low but two peaks were evident energies 
5.7 and 7.5 Mev, shown Fig. 


POSSIBLE REACTIONS 


natural neon contains three isotopes the proportions 90.51% 
0.28% 9.21% and since the induced transitions may excited 
levels the product nuclei, from the energetic point view there are numerous 


TABLE 
(Light) 
(Light) particle 
Energy particle rangein Number Cross 
Level in residual release, energy, chamber, of section, 
Reaction nucleus Mev Mev cm events Loss millibarns 
(A) Type Ne(y, a)O'*, Q = —4.75 Mev 
17.63 Mev Ground level 12.88 10.3 2.45 200 200 
2 6.06 and 6.14 6.8 5.4 .82 3,130 .067 .67 
3 6.91 and 7.12 5.9 4.7 -65 8,550 056 1.80 
4 9.58 (broad) 3.3 2.6 .29 
5 9.84 3.0 2.4 .26 
6 10.36 2.5 2.0 .23 
7 11.25 (broad) 1.6 1.3 .13 
14.8 Mev Ground level 10.0 8.0 1.55 230 -086 
12 6.06 and 6.14 3.9 3.1 .37 2,800 .030 About 1 
13 6.91 and 7.12 2.9 2.3 25 
(B) Type Ne?%(y, a)O'4, = —9.66 Mev 
1 Ey = 17.63 Mev Ground level 8.0 6.4 1.08 360 0.89 .76 
2 1.98 and 6.0 4.8 -69 
perhaps others 
14.8 Mev Ground level 5.1 4.1 
12 1.98 3.1 2.5 28 
(C) Type —12.87 Mev 
17.63 Mev Ground and low 4.76to 4.5 5.93 61,300 
levels at .110, 4.56 
.197 Mev 
2 1.35 3.4 3.2 3.28 About .24 About 5 
20,000 
3 1.57 3.2 0 2.94 
4 2.82 1.9 1.8 1.25 
11 E, = 14.8 Mev Ground and low- 1.9 1.8 1.25 
lying levels 
(D) Ne(y, 2a)C!2, Q = —11.9 Mev 
1 Ey = 17.6 Mev Ground level 5.7 
11 Ey, = 14.8 Mev Ground level 2.9 


(E) Ne(y, ap)N'5, Q = —16.9 Mev 
(F) Ne2(y, Q = —16.91 Mev 
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possibilities for the energy which may released the chamber kinetic 
energy charged particles. assumed that the energy lost per ion 
pair produced recoil fragment not much different from 
the same quantity for protons alphas, then the pulse amplitude produced 
photodisintegration the chamber directly proportional the 
kinetic energy release. pronounced difference between the value for 
light particles and the recoil fragments would tend broaden the peaks. 

Table are listed those reactions resulting charged particles which 
were most likely observed this experiment according the current 
energy level diagrams (Ajzenberg and Lauritsen 1955), together with other 
data relevant the gridded ionization chamber results. 


EXPERIMENTAL RESULTS 


(a) Reactions Observed 

After preliminary runs had established the general form the pulse height 
distributions and counting rates expected, the pulse spectrum was 
examined detail three sections with the kicksorter resolution each 
being made compatible with reasonable running time. Figure shows 
complete spectrum above Mev and the corresponding time-dependent back- 
ground, taken with the gridded ionization chamber and low kicksorter 
resolution. Figure shows the pulse distribution after background subtraction, 
due allowance being made for the extra spreading the two contamination 
peaks which would arise when the gamma flux passes through the chamber. 
The assignment the groups fits well with linear energy scale based the 


position the plutonium alpha peak 5.16 Mev. 
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While there little danger ambiguity the assignment the appro- 
priate reaction the higher-energy peaks, there considerable uncertainty 
the low-energy region the spectrum which shown greater detail 
Fig. 


1500 


1000: 
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5 


ENERGY RELEASE Mev 


_Fic. Pulse height spectrum low-energy particles resulting from the photodisintegra- 
tion neon. 


Firstly, the strong peak marked C2, attributed centered 
3.38 Mev corresponding transition the 1.35-Mev level There 
marked spread this peak towards lower energies might expected 
there were transition similar intensity the 1.57-Mev level 
Secondly, the filling between and Fig. would seem arise 
from reaction A12, i.e. the residual nucleus being left 
the 6.06- 6.14-Mev states. Some might also arise from 
the doubtful level 0.9 Mev (Seale 1953) exists. addition the 
distribution shows evidence for other groups 2.93 2.0, 
1.84, and 1.6 Mev (C4 and/or which might attributed reactions 
indicated Table but for which reliable cross section can estimated. 
The yields from reactions and are expected small because 
the low barrier penetrabilities, hence most the counts between and 
Mev are probably due reaction 413, i.e. proceeding 
the levels near Mev. 

Again the large proton peak was examined under high kicksorter reso- 
lution and with reduced gamma-ray flux appeared that this group was 
rather broader than the others. The over-all stability and chamber resolution 
were hardly adequate produce convincing evidence, yet, can judged 
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from Fig. which consists the addition four separate runs, believe 
that this peak multiple one and could resolved into three components 
similar intensity separated about 100 kev. 
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Fic. Detail the 4.7-Mev photoproton peak from the reaction arising 
from 17.6-Mev gamma radiation. 


This would suggest that the reaction proceeds roughly 
equal proportions the ground state and the two low-lying levels 
110 and 197 kev 


(b) Cross-Section Calculations 

The cross sections were estimated from the area under the peaks. the 
low-energy region where the peaks overlap the areas were estimated 
drawing curves the expected width and obviously the values deduced 
are much more approximate. wall correction was applied the observed 
number events numerical integration based the range the alpha 
particle proton released and this shown Table the column, 
the fraction events not counted the peak. 

The gamma flux was estimated from the counts recorded the standard 
brass Geiger counter, which was placed some distance directly behind the 
chamber. This, together with the intensity ratio 1.6/1 
0°, defined the flux each component the monitor. Corrections were then 
applied for absorption the chamber walls (about per single wall thick- 
ness) and for the inverse square law, the distance from the target chamber 
axis being cm, giving mean flux believed correct +10%. 
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Thus, finally, 


mean square distance neon the sensitive volume 


The cross-section values for the higher-energy peaks are considered 
accurate +20% while the relative values should considerably more 
accurate. For the sake consistency Table contains data obtained only 
with the gridded ionization chamber. The smaller proportional counter con- 
taining pure neon lower pressure gave values within 20% the values 
listed, but was useful only for reactions which released more than 
5-Mev energy the counter. 


(c) Variation Cross Sections with Gamma-Ray Energy 

view the relatively large photodisintegration yields, particularly from 
the reaction, obtained with the gridded ionization chamber, 
attempt was made examine the variation cross sections with change 
gamma-ray energy. using thin lithium targets, less than kev the 
vicinity the 440-kev resonance and about thick away from the 
resonance, the energy the sharp gamma-ray line could taken from 17.59 
18.03 Mev the proton-bombarding energy was raised from 400 900 
kev. The ratio varies quite sharply the 440-kev resonance 
(Devons and Hine 1949) but unfortunately the behavior not known very 
precisely that the interpretation the monitor counter data was some- 
what uncertain. higher proton energies towards Mev there also 
considerable yield 4-Mev radiation from (Kraus 1954), but 
the thick-walled brass counter relatively insensitive this. 

Three parts the photodisintegration pulse spectrum were examined 
corresponding the 4.7-Mev group, Cl, arising from the 
6-Mev group, arising from and the region between 2.5 
and Mev attributed reaction excited state 
transition). The photodisintegration yields per microcoulomb bombarding 
protons the vicinity the resonance for the and groups followed 
moderately well that expected from the intensity the 17.6-Mev radiation 
component according Devons and Hine (1949). sharp fluctuation 
these yields was observed the radiation energy was raised Mev, 
indicating that sharp photodisintegration resonance occurs this energy 
range. However, the data showed that the relative size the group 
fell quite significantly the gamma-ray energy was raised. 

photodisintegrations were observed when the chamber was irradiated 
with and 7-Mev gamma radiation resulting from the proton bombardment 
fluorine targets accord with the data Woods (1952), who showed 
that o(y, was less than cm? this energy. 


DISCUSSION RESULTS 
The alpha spectrum shows evidence for any alpha groups which might 
correspond transitions low-lying excited states confirming the 


t 
t 
i 
+ 
t 
| 
f 
f 
E 
§ 
| 
t 


1164 CANADIAN JOURNAL PHYSICS. VOL. 37, 1959 


current view that the lowest excited level this nucleus lies 6.06 Mev 
above the ground state. 

The values found for the (y, cross sections are agreement with what 
might expected according the selection rules deduced Gell-Mann and 
Telegdi (1953), the basis the charge invariance nuclear forces, for 
which the ground state state and all levels excitation 
about Mev have isotopic spin zero. Electric dipole absorption radiation 
the energies used this experiment would result levels with 
Direct transitions from such levels alpha emission the ground state 
are forbidden the isotopic spin selection rule. Absorption into 
cannot proceed the ground state nor the first-excited state 
6.06 Mev which also 0*, via alpha emission, without violating 
the conservation angular momentum and parity. Thus ground 
state transition these energies proceeds most likely absorption through 
compound state with The formula Weisskopf (Blatt and 
Weisskopf 1952), which provides estimate the upper limit expected 
for the cross section for absorption radiation 17.6-Mev energy, gives 
1.3 millibarns. The experimental value 0.05 satisfactorily less than this 
theoretical limit. also noted that the cross section for 
considerably greater, would expected since isotopic spin 
selection rules impose restrictions this case. 

Again transitions the excited levels 6.14 Mev 
6.91 Mev and 7.12 Mev which are all levels, cannot 
spin conserved. However, these levels are accessible from states 
formed and absorption respectively. The Weisskopf 
formula gives, for absorption, cross-section limit 8.5 millibarns for 
17.6-Mev radiation. The larger cross section observed compared with the 
ground state transition presumably the result the absorption process. 
The observed factor three larger intensity the transitions the 7-Mev 
state compared with the state may result the smaller 
orbital angular momentum the emitted alphas the transitions the 
7-Mev states. very unlikely that significant fraction the counts 
the peak corresponding energy release 5.9 Mev the chamber arise 
from reaction D1, i.e. this forbidden the same way 
the transition the ground state 

The mean cross sections for the broad 14.8-Mev component the 
radiation show the same trends and are similar magnitudes those for 
the sharp 17.6-Mev component. The actual cross sections are little larger 
for the 14.8-Mev radiation, indication perhaps the general trend 
(y, cross section with gamma-ray energy this region. 

The large cross section must arise from transitions from 
levels following absorption. The observed cross-section value 
would suggest cross section for the inverse reaction, with about 5-Mev 
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proton-bombarding energy, 0.4 millibarn, which the order magnitude 
that might expected from reasonable extrapolation Sinclair’s (1954) 
data. There presumably comparable cross section for the 17.6-Mev 
radiation, although the for the neutron reaction —16.9 Mev. would 
expected too that the process would the main method 
breakup the compound nucleus, rather than the (y, process 
observed this experiment. 

These results provide rather satisfactory confirmation the operation 
isotopic spin selection rules the photodisintegration The 
cross section about times large the cross section, result 
similar that found the photodisintegration other light nuclei. 
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ELECTRICAL BREAKDOWN XENON AND KRYPTON 
ULTRAHIGH 


AND 


ABSTRACT 
Breakdown electric fields have been measured resonant cavity for xenon 
and krypton gases the pressure range 0.05 mercury, 
frequency 2800 Mc/sec. Extensive precautions were taken ensure gas 
purity. 


INTRODUCTION 
Measurements the electric fields which cause breakdown gases have 
been used during the past decade investigate the nature the collision 
processes which occur when high-frequency field applied gas (Herlin 
and Brown 1948; MacDonald and Brown 1949; MacDonald and Betts 1952). 
This paper presents experimentally determined values breakdown electric 
fields for xenon and krypton gases high purity. 


EXPERIMENT 


The experimental apparatus similar that used previous investigations 
(Herlin and Brown 1948; MacDonald and Betts 1952). block diagram the 
microwave equipment shown Fig. Power generated magnetron 
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Fic. Block diagram the microwave apparatus. 


operating 2800 Mc/sec travels through coaxial transmission line 
cylindrical cavity resonant the mode. The power divider permits 
variation the power reaching the cavity, and the directional coupler diverts 
small, known fraction this power thermistor and bridge. The direc- 
tional coupler, thermistor, and bridge are replaced slotted section and 
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microwave receiver when the cavity being measured. During the actual 
breakdown measurements, the attenuator, crystal, and meter, which are 
calibrated against the thermistor immediately before each run, are used 
indicate power levels. The electric field strength determined from measure- 
ments the power incident the cavity, the cavity and the voltage 
standing wave ratio resonance, along with the known field configurations 
(Brown and Rose 1952; Rose and Brown 1952). 

The resonant cavity which breakdown takes place must highly 
evacuated before admitting the gas, order ensure gas purity. The vacuum 
system used represented block diagram Fig. incorporates many 


GAUGE 


CAVITY 


DIFFUSION 
PUMP METAL 


VALVE 


FOREPUMP 


COPPER TRAP 


AUTOVAC GAUGE 


COPPER TRAP 


Fic. Block diagram the vacuum system. The dashed line encloses those parts which 
are baked the oven. 


the techniques suggested Alpert (1953, 1954). The resonant cavity 
made oxygen-free high-conductivity copper and connected the pyrex 
system through Kovar seals. Isolation from the pumps accomplished through 
the use metal valve, connected seals and capable 
withstanding baking about 400° Copper foil traps are used instead 
cold traps. The ionization gauge the type, and ion pumping 
used obtain the lowest pressures. The pressure the gas during the 
breakdown measurement determined means Autovac hot wire 
gauge. This gauge used because does not introduce vapors into the 
system would manometers McLeod gauges. Calibration the hot wire 
filaments against McLeod gauge, using sample the gas under study, 
was performed before and after each breakdown run, order take account 
the changes occurring the filaments during the baking the vacuum 
system. 

The pressure the system, before introduction the gas, could kept 
less than mercury, while isolated from the pumps. The gases 
used were reagent grade, supplied the Air Reduction Company. 
sample from each flask used was taken during the breakdown measurements 
and analyzed mass spectrometer (Consolidated Electrodynamics Model 


21-620). 
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The only impurities found the samples krypton used the measure- 
ments labelled Fig. were quantities near the lower 
limit detectability the mass spectrometer. These impurities were: 


less than parts 10°, 
than parts 10* (in one sample only). 
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Fic. breakdown electric field for krypton gas frequency 2800 Mc/sec and 
characteristic diffusion length, 0.101 cm. 

Fic. The breakdown electric field for xenon gas frequency 2800 Mc/sec and 
characteristic diffusion length, 0.101 cm. 


Impurities the gases used for the data Fig. were 
again most cases near the limit detectability, but there was more varia- 
tion from sample sample. Impurities found were follows: 

and addition, 
one sample contained traces and Kr; 
third sample contained less than and less than 
10° 


other impurities were found. 
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RESULTS AND DISCUSSION 


The breakdown field, that is, the electric field which makes gas start 
glow and become conducting, presented Fig. for krypton and Fig. 
for xenon. both cases the field shown for pressures from 0.05 100 
mercury. The cylindrical cavities which the experiments were done had 
height and characteristic diffusion length, 0.101 cm. 
The measurements were made 2800 Mc/sec, which was the resonant 
frequency the lowest magnetic mode the cavities used. 

both figures the solid lines represent the mean values several sets 
data samples the high-purity gas. For these measurements preliminary 
attained before introducing the gas the system. The dashed lines 
both figures represent the data for sets measurements made system 
which was only partly bakable and which used alcohol dry ice vapor traps. 
(Liquid nitrogen traps cannot used because the high boiling temperatures 
xenon and krypton.) The preliminary evacuation this case was 
pressure about mercury, and the pressures the gas during 
the breakdown runs were measured with Octoil manometers. not have 
mass spectrometer analyses for this latter case, but may assumed that 
there was some contamination the gas the system, because the 
conditions described, even though reagent grade gases were used. 

noted that the impurities appreciably alter the breakdown field 
high pressure the case xenon, and lesser extent for krypton, 
while pressures lower than about mercury there effect. 
This difference the breakdown fields may caused air impurity, but 
consideration the ionization potentials and known reactions the major 
components air does not suggest process whereby might happen. 
any rate the experiments show that obtain the breakdown fields character- 
istic pure krypton and pure xenon, the system which the measurements 
are made must completely bakable and not contain sources vapors. 

The electric fields have been measured accuracy better than 
all cases. noted above, was necessary use Pirani-type gauge 
the measurement pressure when obtaining the high-purity data. both 
the high-pressure and the low-pressure ends the range this gauge becomes 
inaccurate. Between 0.5 and 5.0 mercury, pressures were measured 
for each point. Below this range each point has accuracy about 
10°% while above each point has accuracy about 20%. Because the 
curves drawn represent the averages several sets points, the over-all 
accuracy somewhat better and estimate that the error the breakdown 
curves does not exceed for pressures less than mercury, nor 10% 


for pressures above this value. 
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THE 0-0 AND 1-0 BANDS THE SYSTEM NH! 


ABSTRACT 


molecules have been produced during the flash photolysis HNCO, and 
the and 1-0 bands the system have been photographed 
absorption the third order 21-foot concave grating spectrograph. The 
intensity distribution the bands showed that the molecules were approximately 
room temperature. The lines lowest value are identified for branches 
the band and branches the band. The analysis the band 
leads correction the published line assignments for this band. The spin 
splittings both the and electronic states, and the A-doubling the 
state, are considered detail. Wavelengths air are given for lines astro- 
physical interest. 


INTRODUCTION 


The 3360-A bands, the molecule were first described 
Eder (1893), and have been the subject many subsequent investiga- 
tions. Funke (1935, 1936) has published extensive assignments for the lines 
the and bands. Recently Pannetier and Guenebaut (1958) have 
observed the and bands the system emission from atomic 
but with only moderate resolution. 

The determination accurate vibrational and rotational constants for this 
system unfortunately hampered for several reasons: (i) The internuclear 
distances and vibration frequencies the two electronic states are very 
similar, and therefore bands with are very weak account the 
principle; addition the branches the and 
bands are confined very narrow spectral regions with consequent over- 
lapping, and difficult make unambiguous assignments. (ii) When the 
spectrum observed emission from flame very difficult observe 
this wavelength. (iii) Although the branch assignments the bands 
can checked satellite branches (AV are observed, the satellite 
branches are very weak except the lowest values, account the 
rapid transition the state from Hund’s coupling case (a) case 
high-temperature spectra the low levels are not highly populated, and 
very few satellite lines have been observed. Difficulties (ii) and (iii) could 
overcome the bands were obtained with low rotational tempera- 
ture. 

the present work molecules have been produced room tempera- 
ture and the and 1-0 bands have been photographed absorption. 
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EXPERIMENTAL 


molecules were produced during the flash photolysis isocyanic acid 
vapor, HNCO, using apparatus similar that described Callomon 
and Ramsay (1957). The spectrum was observed under the same experimental 
conditions the NCO spectrum (Dixon 1959). The absorption produced 
6-meter path length through the reaction tube was photographed the 
third order 21-foot concave grating spectrograph using Eastman Kodak 
plates, with resolving power about 200,000. The 0-0 band the 
system was extremely strong, and the band appeared 
with about the intensity the band. The 1-0 band was partially 
overlapped weaker band due the NCO radical, but the lines the 
two bands were easily identified account the large ratio (~40) the 
rotational constants the two molecules. few lines the 
component the 3064-A band were also observed, but these lay outside 
the region the band. 

hollow-cathode iron arc lamp was used for wavelength calibration. The 
wavelengths the reference lines were taken from the M.I.T. tables (Harrison 
1939) and used with Edlén’s (1952) vacuum corrections. One plate the 
band and two the 1-0 band were each measured twice comparator 
equipped with photoelectric scanning device similar that described 
Tomkins and Fred (1951). The relative accuracy the frequencies strong 
sharp lines +0.03 


ROTATIONAL ANALYSIS 


Theory predicts that band should have total branches, 
and their lines lowest value are illustrated Fig. these, nine are 
main branches (AV and correspond the nine strong 
branches allowed for band. For the state the 
transition from Hund’s coupling case case very rapid with in- 
creasing rotation, A/B small —2); the remaining satellite 
branches are therefore very weak except the lowest values. previous 
analyses only two satellite branches have been observed (Funke 1935, 1936). 
the present investigation the first line has been identified for the 
branches the band, and for the branches the 1-0 band. 

The lines medium value the and branches the band 
form close triplets and were easily identified. With the values ob- 
tained from these branches was possible assign the and branches, 
and then the and lines low The assignments are indicated Fig. 
The branches the band are more degraded than those the 0-0 
band and were readily identified with the assistance satellite lines, even 
though branches were observed. The assignments for the two 
bands are given Tables and II. For the band these differ consider- 
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ably from the assignments given Funke (1935, 1936), notably for the 
branches and for the and branches low value. the basis the 
present assignment was possible reassign Funke’s measurements the 
high lines. These are given Table III, and will used below the 
discussion A-doubling. comparison the rotational structures the 
band reported Funke and the band reported here shows that 
Funke has correctly analyzed the band apart from few errors low .V. 
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Fic. level diagram and first lines transition. 
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Fic. The 0-0 band photographed absorption the third 
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CONSTANTS FROM THE 0-0 BAND 
molecule state has three levels for each value, with V+1, 
theoretically Schlapp (1937) and the complete expressions for the rotational 
term values are given Herzberg (1950) as: 


(1) 
(2) 

(3) 
where 


and are effective rotational constants, and and describe the spin 
splitting. Thus, the levels follow the simple rotator formula. The 
constants and were therefore determined the usual 
manner, using least squares, from the band combination differences 
obtained from the relations: 


The spin splittings and were each obtained from six different 
relationships between the branches the band, either the separation 
between two adjacent lines the difference between two combination 
differences (see Fig. 1). From equations (4) and (5) seen that: 


and experimentally this sum was found essentially constant. and 
were then determined using iterative process. The calculated splittings 
for each value are compared Fig. with the mean the experimental 
values obtained from unblended lines. 

The lower state combination differences for the band agreed with 
for the band within the experimental error. The rotational constants for 
the level are collected Table IV. 


TABLE 
Constants for the state 


the reorientation the electron spin (AV +1) give rise 
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SPIN SPLITTING 


Fic. observed spin splitting the level the state NH. 


microwave absorption spectrum (Strandberg, Meng, and Ingersoll 1949). The 
transitions are made possible the magnetic dipole moment molecule 


triplet state. should have similar spectrum, and the frequencies 
the lines would given and Figure shows that the 
strongest lines would around (30,000 Mc/s) the rotational 
temperature the molecules were not too high. 


The state tends toward Hund’s coupling case low 
and inverted. Thus there are three substates and the levels 
which can designated and respectively. Each these levels 
exhibits A-doubling. For the labelling the A-doublets shall follow the 
and nomenclature proposed Mulliken (1931) and indicated Fig. 

(1936) has given approximate expressions for the rotational term 
values state intermediate coupling, but these were found 
inadequate for the representation the levels moderate value, 
where the rotational distortion becomes significant while the state has not 
yet approached case Gilbert (1936) has considered the state 
including the effects rotational distortion, but neglecting A-doubling, and 
has shown that the term values the three levels for each are the roots 


third-order secular equation: 
(8) 


Gilbert gives expressions for the elements terms and the rotational 
constants B,, D,, and A,. the unitary transformation matrix which 


diagonalizes 
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The Determination and 
Although not possible give simple expressions for the individual 


term values Gilbert has shown that their sum obeys: 


and therefore determined from the sum the three combination 
differences using the relationship: 


where the are the mean values and 

Many the lines the band are overlapped and the following pro- 
cedure was used determine the constants the levels. The term values 
the upper state were determined adding the term value the appro- 
priate lower level, calculated from the ground state rotational constants 
Table IV, the measured frequency each line the band. These were 


TABLE 
Constants for the state 


(a) Constants using the mean energies the A-doublets 


Constants the A-doubling, and revised values and when 
A-doubling effects are included 


v=0 


averaged for each level, using only those values determined from unblended 
lines. The upper level combination differences were then formed from these 
term values, and the and Table V(a) were determined the method 
least squares using equation (10). 

The Band Origin 

The band origin was derived from the upper level term values summing 
over all levels with the same and using the above determined values for 
the assumption that the A-doubling does not influence equation (9). 

The Spin-Coupling Constant 

The coupling constant for the interaction between the spin and electronic 
angular momenta state may determined for each value 
from equation derived Gilbert: 
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was calculated for each value from the mean energy the and 
levels. The values obtained for were not constant with 
seen Fig. 
35.0 
WITH Y'=0.04 


~~ GILBERT 


50 


Fic. _The spin coupling constant for the level the state, showing 
the variation with when terms containing are neglected. 


Mulliken (1930) has suggested that addition the spin splitting 
state given the coupling constant there can further spin-splitting 
terms follows: 


the upper level term values are modified take account such terms 
found that value that essentially constant with can obtained 
with (Fig. 4). However, this not necessarily the correct 
value for the level, calculated from the mean energies the 
A-doublets. This will discussed below. 


A-Doubling 

The finer interactions state with other states the molecule con- 
tribute extra terms and added the rotational energy. 
general these cause splitting the degeneracy the and levels, giving 
A-doubling: 


34.5 
33.5 
(15) 
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For the low values these doublings were derived directly from the observed 
upper level term values, and high from measurements 
Table III the relationship: 


The interaction *II state with state has been considered detail 
Hebb (1936). The A-doublings each substate can 
given terms three constants and and the elements the 
transformation matrix equation (8). Here, the suffix refers the 
suffix refers the elements appropriate the term value When 
Hebb’s equations were applied was found that for there 
was need higher order term and extra term was 
accordingly added. For his calculations Hebb has grouped together levels 
with the same symmetry for each value, and for his equations must 
use: 


Thus the A-doublings are given by: 


The constants Ci, and were determined using equation (20) the 
method least squares, making use the summation: 


Fig. the experimental A-doublings are compared with values calculated 


20.0 
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according equation (20) with the constants Table The agreement 
good, taking into account that Funke’s measurements some the 
branches merged high values. 

the and the states, which are separated have 
the same quantum number total orbital angular momentum, and are 
related the manner that Van Vleck (1929) describes precession”, 
then the A-doubling constants are given by: 


(22) 
(23) 

(24) 

where: 

(25) 
(26) 
(27) 


interaction constant which Hebb claims always positive, 
and contribution the A-doubling arising from the interaction 
with and states. From the experimental values and and equa- 
tions (23) and (24) obtain —0.132 and 0.0318 com- 
pared with and 0.0358 calculated from equations 
(26) and (27) with The relationship “‘pure therefore 
approximately true. The calculated value extremely small, 
and therefore almost entirely given 

The Spin-Splitting Constant 

the calculation using equation (11) was found necessary 
value that was invariant with the levels lie approximately 
half-way between the and levels the same their term values have 
little effect the calculated value order investigate the spin- 
splitting terms equations more fully, with particular refer- 
equation (8) with the above determined values Dj, and The differ- 
seen that neither the curve for nor that for have the expected 
dependence although the separation between the curves for and 
does increase with Now the the level 
not only will this have marked effect the residuals Fig. for 
but will also lead incorrect value therefore affecting the calcu- 
lated values Fi(.V) and From Hebb’s (1936) equations can 
derive the contribution arising from the interaction 
the state with states. addition, Van Vleck (1932) has shown 
that for Hund’s case (a) the component can interact with and 


— 

— 
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Fic. The difference between the mean energy the observed rotational term values 

and for the level NH, and the values calculated: accord- 
ing equation (8), the effect equation (28) with The full 


states, giving the contribution (equation 22). These same states also 
the component interacts with only states, with only 
for the mean energy and 


the present case can neglected, and the last two terms 
contribute the effective and values, and can ignored. Further, 
the term and the term with coefficient contribute maximum 
0.02 Thus the term involving dominant. The introduction 
predicted equations 14. was found that with —1.6+0.3 
and this could achieved simultaneously for all 
three levels F;(.V), and that agreement with the value 
postulated above (see Fig. 

the main contribution Cy, then and are both negative. One 
might infer from this that the state lies above the 
state, but this conclusion far from being unambiguous. 

Thus possible derive expressions that will represent the observed 
mental accuracy. These involve the nine rotational constants 
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Co, Ci, C3, and xo. Although this sufficient for 12, Funke (1935) 
has shown that for additional term needed, where 


CONSTANTS FROM THE 1-0 BAND 


The Lower State: and 

The lower state vibrational interval and rotational constant 
were obtained plotting the difference frequency between corresponding 
lines the band (these measurements) and the band (Funke 1935, 
1936) against Differences from the various branches agreed 
within the experimental error, indicating that there significant difference 
between the spin splitting the two levels and Values the 
difference were obtained for and within this range the plot gave 
straight line, indicating that Dj’. From the intercept the line, 


The Upper State 

and 

The differences frequency between lines the same the and 
bands were found vary from one branch another. Thus the spin 
splittings the and levels are not identical. The upper level 
term values for were therefore determined the manner described 
above for the level, and the following relationship was used (cf. 
equation 9): 


plot the experimental values the L.H.S. this equation for 
against J(J+1) had definite curvature; therefore Dj. The 
determined the method least squares. 

The Spin-Coupling Constant 

The constant was derived from the experimental term values the 
manner described above for using equation (11). was again found 
necessary postulate spin-splitting constant +0.04 order 

A-Doubling 

The A-doubling the level only available for from the 
present analysis the band. However, values for were obtained 
from Funke’s measurements the band, using equation (16). These 
two sets values agreed the region overlap, and showed dependence 
very similar that for the level. The constants C2, and 
were determined the methods outlined above and are given Table 
seen that the constant C2, which determines the main contribution 
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the doubling high appreciably lower than for the level. 
the same order magnitude (see equations and 27). 


and 

The differences between the observed term values for (mean 
and levels) and those calculated from equation (8) with the constants 
determined above (Table V(a)) showed the same dependence 
illustrated Fig. for Data were only available for and 
within this range the same constants and —1.6 
determined for the level could explain the deviations for 


APPLICATION ASTROPHYSICS 


Should molecules exist interstellar space they could detected 
their ultraviolet absorption. These molecules would have rotational tem- 
perature close and the majority would therefore the lowest 
rotational and vibrational level the ground state (Herzberg 1950, 496). 
thus interest tabulate the wavelengths air those transitions 
that can arise from this This done Table for the 


TABLE 


Wavelengths air and rotational line strengths for the lines the 
system arising from the lowest rotational level 


Wavelength air, 
Rotational line 
Line 0-0 band 1-0 band strength, 


11.6 
2.2 


*The relative vibrational transition probability of the 1-0 band is of the order of 
2% of that of the 0-0 band. 


and bands the system. Not all these lines have 
been observed directly and free from overlapping the laboratory, but their 
wavelengths could calculated using other measured line positions and mean 
lower level combination differences. The relative rotational line strengths were 
calculated according the approximate equations given Budé (1937). 
direct measurement was made the intensity ratio the and 0-0 
bands; this was estimated from the rotational energies the last lines observed 
the and branches the two bands the assumption that the rotational 
distribution approximated that for room temperature. The value 
Table therefore only order magnitude. 
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THEORY THE 
PRESSURE-INDUCED ROTATIONAL SPECTRUM 


ABSTRACT 


The theory induced infrared absorption developed previously applied 
the pressure-induced rotational spectrum hydrogen. The intensity the 
rotational band due mainly the quadrupolar induction effect, and small 
interference effect between the quadrupolar and overlap moments. From the 
experimental data the binary absorption coefficients, values the angle- 
dependent overlap moments are obtained for 
and calculation the overlap moment for pure presented. Rosen- 
type wave functions appear inadequate for calculation the small 
angle-dependent rotational well vibrational overlap moments. The tempera- 
ture dependence the binary absorption coefficient calculated, taking into 
account the quantum effects the pair distribution function, and found 
good agreement with the experimental data. The dependence the ortho— 
para ratio also discussed. The double rotational line has been 
observed and its intensity measured. 


INTRODUCTION 


The pressure-induced rotational spectrum hydrogen has recently been 
studied experimentally compressed gases (Kiss, Gush, and Welsh 1959), 
function the density and the temperature the gas, and for various 
para ratios. This spectrum has been investigated independently Colpa and 
Ketelaar 6). Calculations the binary absorption coefficients the 
rotational spectrum, which were based the quadrupolar induction effect 
and classical statistics for the pair distribution function, have been pub- 
lished (Kiss al. 1959; Colpa al. 19586). the present paper the calcu- 
lation the binary absorption coefficient extended lower temperatures 
taking into account the quantum effects the pair distribution function, 
and the contribution the overlap induction effect the intensity dis- 
cussed. 

the present calculations the intensity the rotational spectrum 
assume that the intermolecular forces are independent the orientations 
the molecules. For such purely central forces there coupling between the 
rotational and translational motions the molecules, and the total energy 
the gas equal the sum the rotational energy and the translational 
energy where and are the quantum numbers specifying 
the rotational and translational motions the molecules. According the 
general theory induced absorption (Van Kranendonk 1957), the spectrum 
then built additively out the separate rotational lines which arise 
from definite rotational transitions accompanied arbitrary changes the 
translational motion the molecules. The true rotational lines correspond 
transitions which the rotational quantum number least one 
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molecule changes such that the total rotational energy different before and 
after the absorption E,). Transitions which only the translational 
energy changes and which the total rotational energy remains unchanged 
can also take place virtue the strong dependence the induced dipole 
moments the intermolecular separations. These transitions give rise the 
translational branch the spectrum, which has been observed mixtures 
the rare gases (Kiss and Welsh 1959a), and for which evidence has also 
been obtained the rotational spectra hydrogen and hydrogen foreign 
gas mixtures (Kiss 1959). The translational absorption arises from 
pure translational transitions 0), and from transitions 
which the translational jump accompanied orientational transition 
(Am, one more molecules, two equal but oppositely directed 
rotational transitions (AJ; The absorption arising from 
all these transitions will referred the translational absorption. 

The integrated intensity, the pure gas, definite rotational line 
given (Van Kranendonk 1957) 


the total volume the gas; and are the normalized Boltzmann 
factors for the rotational and translational states respectively; 
the total induced dipole moment which function the orientations 
cules. The term P,P, (1) corresponds absorption, and the term 
stimulated emission. The sum (1) runs over all pairs translational 
states for which the rotational line does 
not extend zero frequency, the sum over and (1) may extended 
over pairs translational states, since the terms for which 
then give appreciable contribution this sum. virtue the 
completeness the translational states, one then obtains the following 


expression for the integrated intensity (1) (Van Kranendonk 1957): 


The distribution function defined 


and M,, According the result (2), the integrated intensity 
mass the molecules, averaged over all configurations with the help 
the configurational distribution function (3). 


For the integrated translational absorption theorem analogous (2) 
holds, directly evident from the fact that the intensity for frozen 
configuration would zero this case. The restriction now implying 
essential because the translational branch extends zero fre- 


4 
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quency. The calculation the integrated translational intensity therefore 
more complicated. restrict ourselves the present paper the true 
rotational spectrum and discuss the calculation the translational intensity 


separate paper. 


GENERAL EXPRESSIONS FOR THE BINARY ABSORPTION 
COEFFICIENTS THE ROTATIONAL LINES 


The expression (2) for the absorption coefficient can developed 


powers the density N/V the gas much the same way was 


done previously for the vibrational spectrum (Van Kranendonk 1957). The 
only difference that the vibrational case the vibrating molecule plays 
particular role, whereas the rotational case all the molecules should 
treated par. One must then introduce U-functions means the 


equations 


(12) U(12); 


etc. 


Proceeding this way, the expression for the binary absorption coefficient 


the expansion 


(5) a= ayn? + aon? + 


the absorption coefficient powers the density found given 


rotational line and go(R) the low-density limit the pair distribution 
function g(R). The quantity 

(7) M,,-(R) @2; R)|r’) 


the matrix element the pair dipole moment This pair dipole moment 
expanded terms spherical harmonics the orientations the mole- 


cules, 


(8) Bx (1, R) =4r R) YX} (@1) (we). 


which lies along and (6;, specifies the orientation molecule 
relative polar axis. For homonuclear diatomic molecules for 
odd. further note the symmetry property 


which consequence the invariance the Hamiltonian the two 


— 
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interacting molecules under inversion point midway between the 
molecules, and the relations 


which follow from reality requirements and from the transformation properties 
for rotations around the z-axis respectively. 

When the expansion (8) substituted into (6), the resulting expression 
for the binary absorption coefficient can written the form 


(14) To(AiA2) | | R) * g)(R)dR, 


KM IM? 


where the notation the same the vibrational case (Van Kranendonk 
1958). take into account the terms the expansion (8) corresponding 
(2,0), (0, 2), and (2, 2), and neglect the terms corresponding 
because the symmetry property (9). For the single rotational lines corre- 
sponding the transitions J+2, the binary absorption coefficient (12) 
then given 


The factors L(J) and L’, defined 
J=0 
where 


determine the dependence the absorption coefficients the distribution 
the molecules over the rotational states, whereas the quantities and 
determine the temperature dependence the absorption coefficients 
due that the pair distribution function. For the rotational lines 
corresponding the double transitions the binary 
absorption coefficient 


\ 
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Finally, the total binary absorption coefficient the whole rotational band, but 
exclusive the translational branch (cf. Section 1), given 


J=0 J<J’=0 


The calculation the binary absorption coefficients thus reduced the 
calculation the quantities and defined eq. (14), i.e. 
the calculation the coefficients the expansion (8) the induced dipole 
moment. the case the vibrational spectrum (Van Kranendonk 1958), 
assume that the induced dipole moment consists short range overlap 
moment arising from the overlap forces between the molecules, and long 
range moment due the quadrupolar induction effect. The calculation 
these two parts the induced moment presented the next section. 


CALCULATION THE INDUCED ROTATIONAL DIPOLE MOMENT 


The main contribution the induced rotational dipole moment comes 
from the quadrupolar induction effect (Van Kranendonk 1958). The intensity 
the rotational spectrum resulting from this induction effect has been 
calculated Colpa al. (19586) and Kiss a/. (1959). our notation, 
the non-vanishing coefficients the expansion (8) for this part the induced 
moment are given 


(23) 


plus other coefficients which follow from (23), (24) means the 
symmetry relations (9) and (10). the quadrupole moment, 
the average polarizability, and the anisotropy the polariza- 
bility hydrogen molecule the normal state. The phases the spherical 
harmonics (8) have been chosen rather than The spherical 
The rotational overlap dipole moment for pair hydrogen molecules 


given 


(1234) the electron wave function function the position 
and spin co-ordinates the four electrons, depending parametrically the 
nuclear co-ordinates. The sum over (24) runs over the four protons 


= 
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and the position vector the ith particle. For use the same type 
wave function was used for the vibrational case, viz. 


where the operator exchanging the electrons and x(12) the singlet 
spin function, and 


(29) 


The one-electron orbitals are distorted Rosen-type wave functions (Van 
Kranendonk and Bird 1951), 


(31) = Kant act 


etc. The term containing describes the distortion orbital due atom 
the terms containing and describe the distortions induced the 
atoms and respectively. The orbitals and are defined similarly, 
with 


(32) | 1, = = Kea—Aac— Ape; 


The dipole moment (25) vanishes the overlap between the two molecules 
and the distortions described the terms containing are neglected. 
calculate the moment (25) the first non-vanishing order the square 
the overlap and the X’s, treating the squares the overlap integrals and 
the parameters small quantities the first order. The resulting expression 
for equal the sum the following two terms: 


(33) Uaet 14S 


The summation signs indicate that three additional terms should written 
corresponding cyclic permutations and The quantities 


Sac 


appearing (33) and (34) are evaluated putting all equal 
zero, whereas the quantities and (33) contain only the terms 
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and respectively that are linear the \’s. The moment (33) there- 
fore proportional the and called the deformation moment, whereas 
the moment (34) proportional the square the overlap integral between 
the two molecules, and called the exchange moment. The quantities (35) 
can easily evaluated, giving 


(36) a (2? / drs; 


and similarly for etc. The moments (33) and (34) have been evaluated 
molecules, and for number intermolecular separations The resulting 
moments turn out decrease very nearly exponentially with increasing 
and lead the following values the expansion coefficients (8): 


where 2.93A the molecular diameter appearing the 
potential. The coefficients (37) should added the corresponding co- 
efficients (23) arising from the quadrupolar induction effect get the total 


expansion coefficients. This completes the evaluation the induced rotational 
dipole moment. 


CALCULATION THE BINARY ABSORPTION COEFFICIENTS 
THE ROTATIONAL LINES 


When substitute the expressions (23), (24), and (37) for the expansion 
coefficients into the expression (14) for the integrals obtain the following 
expressions for the binary absorption coefficients: 

For the single rotational lines, cf. eq. (15), 


For the double rotational lines, cf. eq. (19), 


For the total binary absorption coefficient, cf. eq. (21), 


The dimensionless quantities and are defined (cf. Van Kranen- 
donk 1958) 


Ax ’ 5» (= ’ 


5 
eo eo 


+ 
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and the quantity 
(43) 


has the dimension binary absorption coefficient. The integrals and 
are defined 


0 


0 


‘ Zo( ) ’ 
where high temperatures classical statistics may used, 
and the pair distribution function then equal the Boltzmann factor 
exp where V(x) the intermolecular potential for which 
choose the potential 


(45) V(x) 


intermediate temperatures, go(x) may expanded powers Planck’s 
constant. The resulting expansions the integrals (44) are written the 
form 


(46) 


and similarly for and where The integrals 
and etc. have been evaluated numerically for number reduced 
temperatures The results are collected Table note here 
that Table previous publication (Van Kranendonk 1958) the headings 


TABLE 


Classical values and first and second quantum corrections the radial integrals and 
T* Tot Kei Jel K® yo 
3.0 2.39 0.37 2.23 3.67 
8.0 3.98 0.068 3.12 4.37 
0.0 2.14 2.85 


60 10.76 13.91 —1.37 0.056 


the 4th and the 5th columns this table should and 
respectively. the calculations the correct values and 
were used. 
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first discuss the total binary absorption coefficient 300° Using 
the theoretical value James and Coolidge (1942) for the quadrupole moment, 
mann 1935), the numerical values the various terms the expression (41) 
for the total binary absorption coefficient are given 


where The first term the intensity due the overlap 
moment alone, and quite small. The second and fourth terms are due 
the quadrupolar induction effect, and give the main contribution the total 
intensity. The third term represents the interference effect between the 
quadrupolar and overlap moments. The experimental value given 
Kiss al. (1959) the total intensity inclusive the translational absorption, 
whereas (47) the calculated intensity exclusive the translational absorption 
(cf. Section 1). The experimental value must therefore corrected for this 
This has been done analyzing the shape the rotational absorption 
lines and Welsh and the resulting experimental value 
The theoretical value (47) good agreement with this value 
far the order magnitude concerned, but about 20% too small. The 
quadrupolar induction effect gives contribution 1.06@ and thus accounts 
for 90% the observed intensity. This value depends the quantities 
and all which are known rather accurately for hydrogen. The 
calculated value the quadrupolar intensity should therefore rather 
accurate. The inclusion the overlap induction effect reduces the intensity 
giving final value about 20% smaller than the experimental one. 
possible that this discrepancy due inaccuracy the calculated 
quadrupolar intensity, the fact that use simple model for the total 
induced dipole moment, which only the terms longest and shortest 
(exp range are taken into account. Thus addition the terms, 
there will terms varying etc. is, however, more likely that the 
calculation the overlap moment discussed Section which based 
the use Rosen-type wave functions, not reliable. Support for this con- 
clusion obtained from analysis the vibrational-rotational spectrum 
hydrogen. The most important part the vibrational overlap moment 
independent the orientations the molecules (Van Kranendonk 1958). The 
order magnitude this part the overlap moment given correctly 
the calculation based Rosen-type wave functions. However, the angle- 
dependent part the vibrational overlap moment order magnitude 
smaller than the isotropic part, and not given correctly the calculations. 
detailed analysis the intensities the vibrational spectrum, carried out 
Hunt (1959), has shown that the phenomenological value the angle- 
dependent part the overlap moment, obtained fitting the calculated 
absorption coefficients the experimental ones, the same order 
magnitude the calculated value but opposite sign. the case the 
rotational overlap moment, the isotropic part the moment vanishes identic- 
ally result the symmetry relation (9), and only small angle-dependent 
part survives. When neglect the first term the expression (41) for the 
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total binary absorption coefficient, the contribution the overlap moment 
the total intensity can described means the single parameter 
The expression (41) then becomes 


All the quantities appearing this expression are known, except 
empirical value can obtained adjusting the value (48) the ex- 
perimental value This has been done for 300° giving 
the case the vibrational overlap moment, reach the con- 
clusion that the calculations based Rosen-type wave functions are not 
elaborate enough yield reliable value for the angle-dependent part the 
overlap moment. therefore adopt for the empirical value and 
calculate the temperature dependence and the intensities the separate 
rotational lines, using this value for theoretical determination based 
the use more realistic wave functions for two interacting hydrogen 
molecules, will not attempted the present paper. 

The temperature dependence has been calculated using the expression 
(48) with The quantities and appearing (48) have 
been calculated numerically functions the temperature for equilibrium 
and normal hydrogen, and the results are collected Table II, where and 


TABLE 
Statistical factors and for hydrogen 
Equilibrium Normal 

0.5 1.000 0.500 0.700 0.455 
1.0 0.967 0.500 0.7 0.455 
1.5 0.882 0.493 0.700 0.455 
2.0 0.810 0.481 0.698 0.455 
2.5 0.762 0.470 0.697 0.453 
3.0 0.729 0.460 0.692 0.450 
3.5 0.705 0.452 0.684 0.446 
4.0 0.686 0.444 0.675 0.441 
5.0 0.656 0.429 0.652 0.429 
6.0 0.630 0.415 0.415 
7.0 0.606 0.400 0.606 0.400 
8.0 0.584 0.386 0.584 0.386 
9.0 0.564 0.373 0.564 0.373 
10.0 0.545 0.360 0.545 0.360 


hydrogen well for mixtures hydrogen with foreign gases. The resulting 
temperature variation shown Fig. for the equilibrium 
ratio each temperature. The agreement with the experimental temperature 
variation seen very good. There marked difference between 
the temperature dependence the total intensity the vibrational band and 
that the pure rotational band. For the vibrational band the intensity 
independent the distribution the molecules over the rotational states 
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(Van Kranendonk 1958), and the temperature dependence due entirely 
the temperature dependence the pair distribution function. For the 
pure rotational band this not so; important part the temperature 
variation due the factors and appearing (48). These factors 
depend the distribution the molecules over the rotational states, and 
therefore the temperature. 

The variation the total intensity (48) with varying ratio 
given temperature due entirely the dependence the factors and 
this ratio. Neglecting the small term proportional the ratio 
the intensities for normal and equilibrium hydrogen given 


where given eq. (22). According (49), the intensity liquid air 
temperature should less for normal than for equilibrium hydrogen. 
Experimentally (Kiss 1959) the value for normal hydrogen 85° 
about less, satisfactory agreement with the theoretical value. 

The intensities the separate rotational lines can calculated using the 
expression (39). The calculated intensity ratios can used analyze the 
observed absorption contours; this will discussed detail separate 
paper. liquid air temperature the and S(1) lines are well resolved, 
and the observed ratio the intensities for equilibrium hydrogen then 
equal S(0)/S(1) 1.06, whereas the theoretical value this ratio 1.14. 
The double rotational transition has also been observed, and the 
intensity this line estimated about cm®. The value 
accurate measurement the intensity necessary decide whether the 
difference between the calculated and observed value significant. 

Finally, discuss the rotational spectrum hydrogen induced foreign 
gases (Kiss 1959). The total binary absorption coefficient for mixture 
hydrogen with foreign gas given 


wi 
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assume that the constants characterizing the intermolecular potential 
for two different molecules are given and all 
the quantities appearing (50) are known for the mixtures interest, except 
Values for the various mixtures can determined fitting the 
value (50) the observed values This leads the following 
positive value required. each case the quadrupolar induction effect 
accounts for least 90% the total intensity, and the values given here 
should therefore regarded only order magnitude estimates. This 
true particular for the mixtures, for which the constants and are not 
well known. 

conclusion can say that the theory the induced infrared absorption 
based the (exp —4) model for the induced dipole moment accounts satis- 
factorily for the observed properties the pressure-induced pure rotational 
spectrum hydrogen the pure gas and mixtures. 
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NOTES 


NOTE s-WAVE AVERAGE NEUTRON CROSS SECTIONS* 
KAMAL SETH 


was shown Thomas that the average cross-section problem 
(Feshbach 1955) corresponds the averaging the collision matrix 
over region containing many resonances, and hence consists replacing 
the logarithmic derivative its average The average 
function given 

where related (and sometimes confusingly called) the 
strength function, and (the integral extends over the 
energy interval averaging). For the s-wave neutrons one conventionally 
puts the boundary parameter and cbtains (Lane and 
(putting and omitting the summation over spin states for the sake 
simplicity) 

(1+ 
which naturally the same expression which obtained for the total cross 


section for optical, square-edge, spherically symmetric potential well 
(eq. 3.10 Feshbach al. 1955) the following identity established. 


(1) 


Here and are the real and the imaginary parts respectively the re- 
ciprocal the logarithmic derivative for the potential the surface. 

impractical use eq. the analysis the experimental data for the 
purpose determining and For small values Feshbach, Porter, and 
Weisskopf (1955) approximated eq. 


Lane and Lynn (1957) that this expression neglects the 
interference term which the same order the first term, i.e., the 
order Lane and Lynn (1957) obtained 


expanding both the numerator and the denominator the right-hand side 


eq. and retaining terms orders p?. This expression was used 
preliminary analysis the author (Seth resolve the discrepancy 


*Supported the United States Atomic Energy Commission. 
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between the thin and thick sample results for strength function 
obtained Marshak and Newson (1957) from the analysis the average 
total cross-section experiments for nuclei with atomic weights between and 
75. The energy-independent term the braces eq. (identified with 
Marshak and Newson (1957)) was shown negative for 50, 
and found important whenever the strength function has large value 
(Seth more exhaustive study has revealed that while eq. quite 
accurate (within eq. for energies kev, may gross error 
for higher energies. found that inclusion terms order makes the 
expression unwieldy without improving its accuracy appreciably. order 
able analyze total cross-section data the region Li(p, neutron 
energies simply and with accuracy find the following approximation useful: 


small. This permits further simplification 


cos 


This expression found much more accurate than eq. without being 
any more complicated. Fig. have plotted the errors eqs. and 
expressed percentage the exact cross section eq. The values 
and used are those obtained Feshbach, Porter, and 
Weisskopf (1956) for the Saxon-type optical potential, 
0.45 fermi. kev Lane and Lynn’s (1957) approximation seen 
(Fig. error ~10% for all values reaching maximum 
values 25% and 170. The situation rapidly deteriorates 
with increasing energy. 500 kev the error ~100% almost everywhere, 
becoming high 250% 170. contrast, our approximation, 
eq. seen (Fig. have less than error for all values energy 
the regions and The error other places usually in- 
creases with energy and oscillates between positive and negative values which 
are never greater than +25%, usually within 10%. low energies, kev, 
course eq. equivalent eq. and has the same <1%. Further, 
may noted that not only the magnitude errors our approximation 
smaller but the energy dependence the errors weaker than Lane and 
(1957). Thus Lane and Lynn’s errors increase factor ~35 (from 
170%) over energy interval from kev 500 kev, while, 
over the same energy interval, our errors show almost increase for most 
atomic weights; the worst case they increase factor ~10 (from ~2% 
~20%). 

may mentioned here that the form eq. was suggested the fact 
that simple averaging the single-level formula derived 


x 


NOTES 
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5 j 
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Fic. Percentage error the different approximations eq. for exact s-wave average 
cross section. Lane and (1957) approximation, eq. note the logarithmic scale for 
the ordinate. (b) Our approximation; errors less than +0.1% have not been drawn. 


from the R-matrix theory (with the inclusion the ‘‘background 
whose effect replace everywhere the conventional single- 
level expression) yields exactly the first and the third terms the eq. The 
second term, mentioned earlier, represents the average contribution the 
interference and can not obtained single-level 
approximation. However, its proportionality can easily 
established evaluating the average contribution the corresponding term 
the scattering cross section obtained simply adding two single-level 
scattering amplitudes, suitable averaging more sophisticated two- 
level formula (Lane and Thomas 1958). 


1201 
1000; — 
kev 
2 t 1 
A 
| | 
| | SOO kev/ 
8 4 —— + _ +$——— + + + - 
| J | 
| | | | | 
-24 


1202 CANADIAN JOURNAL PHYSICS. VOL. 37, 1959 


must pointed out that the use any the above expressions 
interpreting data only valid the experiment measures true cross section. 
case the above expressions can only used for analyzing cross sections 
obtained very thin samples where the sample thickness). 
For thicker samples, methods based single-level area considerations have 
been used (Marshak and Newson 1957; Seth 1958a, b). Thomas has 
also developed suitable expression for interpretation average transmissions 
terms the average collision function, and the method, though cumbersome, 
has been used successfully (Darden 1955; Boreli and Darden 1958). 
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NOTES PHOTOGRAPH THE STARK EFFECT 


The accompanying photograph may serve temporarily suitable illustra- 
tion the Stark effect and indeed being used Professor Jenkins 
reference work where the text necessarily limited. The photograph some- 
what enlarged from plate taken with Surdo d-c. source which the 
electric fields accompanying the cathode fall potential are distributed over 
axial distance only millimeters. Owing heavy-particle excitation, 
relatively strong light emitted from the core the otherwise dark space. 

The phenomenon revealed entered physics through the quantum theory 
spectra and the rather intricate web experimental facts interpreted 
quantum mechanics. 

The present illustration shows polarized components the parhelium group 
roughly centered round the diffuse series line \4144. 
Analyses for individual lines are essentially asymmetric states degenerate) 
while the group builds over-all symmetry the external electric field 
increases (mutual repulsion between states with same value.) Incidentally, 
molecular lines the background show little displacements. 
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For the maximum displacements, just under 100 the field 85,000 
v/cm. therefore apparent that few thousand volts per centimeter 
enough cause the necessary displacements and bring out the new electric 
combination lines, corresponding equal any possible value. 

the extreme right member the first series recognized 
having complex pattern for each its members external electric 
fields. The pattern characterized very weak outer component. More 
generally, for any parhelium transition the normal pattern 
possible. 

White dots mark the four, i.e. components corresponding 
transitions from |m| initial states and these components illustrate very 
well the symmetry which develops the displacements become large compared 
with the zero field separations. When odd, however, there 
component and course the state highest |m| value undisplaced 
first approximation. 

the image, appears vanish passes under the zero 
field position Indeed, many components must vanish when the 
displacement reaches the critical value which factor the intensity goes 
zero. 

From text books one may conclude that these well-established features are 
not widely known. Meantime experiments are proceeding under very high 
resolution. 
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LETTERS THE EDITOR 


Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should submitted not later than six weeks 
previous the first day the month issue. proof will sent the authors. 


The Function Constitutional Supercooling Crystal Growth 


This letter comment intended correct one basic misconception presented recent 
paper Kirkaldy growth and the thermodynamics irreversible 
(Kirkaldy 1959). concerns itself only with the conclusions that the author makes with respect 
the existence constitutional supercooling the liquid ahead cellular dendritic 
interface morphologies. Kirkaldy concludes, erroneously, that constitutional supercooling 
eliminated the cell caps necessary consequence the minimum entropy production 
condition. will shown here that such conclusion unwarranted since necessitates 
violation one the boundary conditions the problem. 

The principle minimum rate entropy production non-equilibrium systems, like the 
principle minimum free energy equilibrium systems, presents the optimum thermodynamic 
state the system but does not indicate the probability reaching this state finite time. 
severely segregated freezing crystal exhibits lower rate entropy production than crystal 
freezing homogeneously and therefore preferred (Kirkaldy 1959). However, constitu- 
tional supercooling exists the liquid, the interface remains planar and the crystal homo- 
geneous, the crystal continues freeze non-optimum manner. certain degree 
constitutional supercooling cellular interface forms and segregated crystal produced. 
through the existence constitutional supercooling that the desired end state reached. 
will shown later publication (Bolling and Tiller, published) that cells will only 
form the caps are stable higher temperature than planar interface. Because the 
change equilibrium temperature due curvature, this will only occur when critical degree 
constitutional supercooling reached. 

The interface morphology does not change order eliminate the constitutional super- 
cooling, although the author was once this opinion, the interface morphology changes through 
the aid the constitutional supercooling produce the segregated state lower rate 
entropy production. fact very little constitutional supercooling eliminated when cells 
form. Rather than carry out rigorous treatment the problem may simply shown that, 
for most growth conditions, constitutional supercooling cannot eliminated the cell caps 
without violating the flux condition the cell caps. 

The solute flux condition the cell caps 


where the rate freezing, the solute partition coefficient between the two phases, 
the diffusion coefficient, and the solute concentration the liquid. The concentration 
given where the concentration far from the interface and the 
build-up the interface. For planar interface that for 0.5 
where The condition zero constitutional supercooling the interface given 


where the liquidus slope. Now let and grad adequate for cell formation with 
wt.% Sn. With the same and grad obvious that for equations 
(1) and (2) cannot simultaneously satisfied with wt.%. Thus, constitutional 
supercooling cannot possibly eliminated the caps cells for most growth conditions 
conditions. 

The general conclusion must that cells will adopt that morphology compatible with 
thermal and constitutional constraints that produces the minimum rate entropy production. 
The constitutional supercooling acts bridge for this morphology change. After the interface 
morphology change the constitutional supercooling still remains. 


BOLLING, G. F. and TILLer, W. A. To be published. 
KIRKALDY, J. S. 1959. Can. J. Phys. 37, 739. 
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Zero Spin the 1265-kev State 


known that the first two excited states 328 and 621 kev have spin and 
even parity (Johns and Nablo 1954; Mandeville, Varma, and Saraf 1955). The purpose this 
letter present angular correlation evidence that the third excited state 1265 kev has 
zero spin and, presumably, even parity. The spin sequence, this nuclide therefore 
identical with that (King and Johns 1959) and this zero level represents the second 
example the existence beta vibrations the region intermediate coupling between the 
rare earth nuclei and those near the 126-neutron closed shell. 

Neutron irradiation iridium produces 70-day and 19-hour While most the 
nuclei decay directly the ground state the nuclei all decay via cascading 
gamma rays. Thus even with iridium enriched 89% the coincidences 
from are not negligible. One-milligram samples enriched iridium metal sealed small 
quartz vials were irradiated the Brookhaven reactor for hours and used unopened 
condition. For each cascade interest the angular correlation function was recorded with 
carefully centered sources for nine angles between 90° and 270° using equipment 
described earlier publication (King and Johns 1959). Solid angle corrections, for the 

Nal crystals used, were experimentally determined, using the method de- 
scribed Lawson and Frauenfelder (1953), for series gamma-ray energies from 129 
1330 kev. 

Since both the and 643-621-kev cascades start from the 1265-kev level Pt!™ and 
involve spin sequences the form one might expect the correlation functions for these 
cascades identical ideal experiment. Unfortunately, the resolution scintillation 
spectrometer insufficient make possible clean separation the many radiations present 
and corrections must made for competing cascades. 


2.0 
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Fic, 1, Angular correlation functions for the two cascades depopulating the 1265-kev level in Ir!%. The dotted 
curves show the least-squares fit to the experimental points. The dashed curves are the same corrected for the 
finite solid angle of the detectors. The solid curves are the expected theoretical functions. 


The single-channel spectrum obtained this work was very similar that presented 
Fig. Mandeville’s paper referred above and not shown here. Data was recorded with 
the scintillation spectrometers set register coincidences between the following portions 
the single-channel spectrum: photopeak with itself, 328-kev photopeak with the 
photopeak, and 328-kev photopeak with the entire spectrum above 1000 kev. From 
the results the last-mentioned pair, estimate the contribution the observed 
coincidence rate from competing cascades can made. Figure presents the experimental 
correlation data after the subtraction small correction (<10%) for accidental coincidences, 
together with the least-squares fit the experimental points. 

The figure also shows the experimental curves corrected for the finite solid angles the 
detectors and the theoretical curve expected for 0-2-0 cascade. Before the solid angle 
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corrected curves are compared with theory, they must adjusted for the effects competing 
cascades and The data available shows that the subtraction these effects 
accentuates the anisotropy the curves Fig. Using our best estimates these corrections, 
the value W(130) depressed 0.5 while W(180) raised 1.75. Because the magnitudes 
these corrections are somewhat uncertain, have preferred publish only our experimental 
curves and draw our conclusions from them. 

The two cascades under investigation are both the form where has integral 
value less than The only choice for which gives large positive coefficient zero; none 
the other choices give correlation functions resembling the experimental ones. Thus 
despite the rather poor agreement between the experimentai coefficients and the theoretical 
coefficients for 0-2-0 cascade, one can conclude with confidence that the 1265-kev state has 
zero spin. consideration the intensities the beta rays feeding this level suggests that its 
parity even, though this conclusion cannot drawn from the present experiment. 


These results were presented the fourteenth congress the Canadian Association 
Physicists early June, 1959. The authors would like acknowledge recent communication 
from Dr. Marklund Uppsala describing work which presented the Stockholm 
conference during the same week. His conclusions are agreement with ours. 
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